Method, apparatus, and program for measuring wavelength dispersion of optical waveguide

ABSTRACT

The invention is to provide a wavelength dispersion measuring method for an optical fiber capable of measuring a wavelength dispersion in a section extending from an incident end of the optical fiber to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical fiber is independent of time. Incident light having a known spectral density function S(ω) is caused to be incident on the optical waveguide, and a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by a correlational function between a signal being proportional to a scattered light amplitude generated at a first halfway point and a signal being proportional to a scattered light amplitude generated at a second halfway point.

CROSS-REFERENCE TO THE RELATED APPLICATIONS

This application is a national phase of the international application No. PCT/JP2006/318180 filed Sep. 13, 2006, the entire contents of which are incorporated by reference. This application also claims benefit of priority under 35 U.S.C. §119 to Japanese Patent Application No. 2005-276372 filed Sep. 22, 2005, Japanese Patent Application No. 2006-013796 filed Jan. 23, 2006 and Japanese Patent Application No. 2006-048759 filed Feb. 24, 2006, the entire contents of all of which are incorporated by reference.

TECHNICAL FIELD

The present invention relates to a method, an apparatus, and a program for measuring wavelength dispersion on an optical waveguide.

BACKGROUND FIELD

Since an optical communication system using an optical fiber serving as one of optical waveguides, wavelength dispersion of the optical fiber considerably affects the characteristics of the optical communication system to distort a signal waveform in propagation of an optical signal, it is indispensable in design of the communication system that wavelength dispersion of the optical fiber or the like is known by measurement. In particular, in an actual operation or a test of a system, when an optical cable is installed, distribution of optical fibers in a direction of length of wavelength dispersion must be often measured. In this manner, as a method of non-destructively measuring a distribution of optical fibers in the direction of length of wavelength dispersion, a conventional method of using optical four-wave mixing in optical fibers is given (for example, see Non-patent Document 1).

Non-patent Document 1: L. F. Mollenauer et al. “Method for facile and accurate measurement of optical fiber dispersion maps, “OPTICS LETTERS, VOL. 21, NO. 21, Nov. 1, 1996, pp. 1724-1726

Non-patent Document 2: E. Brinkmeyer and R. Ulrich, “High-resolution OCDR in dispersive waveguides, “Electronics Letters, vol. 26, No. 6, pp. 413-414 (15 Mar. 1990)

DISCLOSURE OF THE INVENTION Problem to be Solved by the Invention

However, since the conventional method must use a nonlinear optical effect in an optical fiber, the method has a disadvantage of requiring a laser beam having a high intensity (in Non-patent Document 1, a measuring system is constructed by a laser beam having about 1 W). For this reason, the present invention provides a new method of non-destructively measuring a distribution of optical waveguides in a direction of length of wavelength dispersion by a method having an essential principle different from that of the conventional method. It is an object to provide a method, an apparatus, and a program for measuring waveguide dispersion without requiring a high-output light source unlike in the conventional technique.

Means for Solving the Problem

The present invention is a wavelength dispersion measuring method for an optical waveguide capable of measuring a wavelength dispersion in a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, including: a step of causing incident light having a known spectral density function S(ω) to be incident on the optical waveguide; a step of observing, by the incident light being incident on the optical waveguide in the step of causing the incident light to be incident on the optical waveguide, a signal being proportional to a scattered light amplitude generated at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and a signal being proportional to a scattered light amplitude generated at a second halfway point being different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way; and a step of calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by a correlational function between the signal being proportional to the scattered light amplitude at the first halfway point and the signal being proportional to the scattered light amplitude at the second halfway point, the signals being observed in the step of observing the signals which are proportional to the scattered light amplitudes.

In the wavelength dispersion measuring method, the step of observing the signals being proportional to the scattered light amplitudes is desirably the step of observing of Rayleigh-scattered lights as the scattered lights at the first halfway point and the second halfway point.

In the wavelength dispersion measuring method, the step of observing the signals being proportional to the scattered light amplitudes is desirably the step of guiding scattered light generated by incident light being incident on the optical waveguide and incident light passing through delay means having delay time equal to propagation time corresponding to twice a length between the incident end of the optical waveguide and the first halfway point to first interference signal detecting means and second interference signal detecting means which have a relative delay difference equal to a distance between the first halfway point and the second halfway point to detect an interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and an interference signal I₂ being proportional to a scattered light amplitude from the second halfway point, and the step of calculating the wavelength dispersion is desirably the step of calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by a correlational function between the interference signal I₁ and the interference signal I₂.

In the wavelength dispersion measuring method, variable delay means set on a path of the scattered light or the incident light is arranged, the step of observing the signals being proportional to the scattered light amplitudes is desirably the step of respectively detecting the plurality of interference signals I₁ and I₂ while causing the variable delay means to change quantity of delay, and the step of calculating the wavelength dispersion is desirably the step of calculating a correlational function between the interference signals I₁ and I₂ by calculating the following numerical expression (40) expressing an ensemble average of products of the interference signals I₁ and I₂ detected in the step of observing the signals being proportional to the scattered light amplitudes: I₁I₂*  [Numerical Expression 40]

In the wavelength dispersion measuring method, variable delay means set on a path of the scattered light or the incident light is arranged, the step of observing the signals being proportional to the scattered light amplitudes is desirably the step of detecting the interference signals I₁ and I₂ as interference signals I while causing the variable delay means to continuously change delay time τ, and the step of calculating the wavelength dispersion is desirably the step of recording the interference signals I detected in the step of observing the signals being proportional to the scattered light amplitudes as a function I(τ) of τ and calculating a correlational function between the interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and the interference signal I₂ being proportional to a scattered light amplitude from the second halfway point by calculating the following numerical expression (42) expressing an ensemble average of the interference signals I₁ and I₂: I₁I₂*  [Numerical Expression 42] by the following numerical expression (41): I ₁ I ₂*=

I(τ_(i))I(τ_(i) +d/c)

_(i)  [Numerical Expression 41] (where < >_(i) expresses an averaging operation about τ_(i)) by using the function I(τ).

In the wavelength dispersion measuring method, by using the following numerical expression (43) calculated as the correlational function of the interference signals I₁ and I₂: I₁I₂*  [Numerical Expression 43] a dispersion D in the section extending from the incident end of the optical waveguide to the arbitrary point along the way is desirably calculated by the following numerical expression (44):

$\begin{matrix} {\mspace{79mu}\begin{matrix} {\frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{1}}^{2}}} = \frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{2}}^{2}}}} \\ {= \frac{\int_{- \infty}^{\infty}{{\Gamma(Z)}{\Gamma^{*}\left( {Z - d} \right)}{\mathbb{d}Z}}}{\int_{- \infty}^{\infty}{{{\Gamma(Z)}}^{2}{\mathbb{d}Z}}}} \end{matrix}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 44} \right\rbrack \\ {\mspace{79mu}{{{where},{{\Gamma(z)} \cong {{\exp\left\lbrack {{\mathbb{i}}\;\frac{\omega_{0}}{v_{p}}z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}}}\mspace{79mu}{D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}}}} & \; \end{matrix}$

-   S(ω): spectral density function of known incident light -   d: relative delay difference between first and second interference     signal detecting means -   ω₀: central frequency of incident light -   ν_(g): group velocity of optical waveguide (can be calculated from     refractive index or the like of optical waveguide) -   z₀: length of section extending from incident end of the optical     waveguide to arbitrary point along the way -   λ: central wavelength of incident light |I₁|²·|I₂|²: squares of     absolute values of I₁ and I₂.

The present invention is a wavelength dispersion measuring apparatus for an optical waveguide capable of measuring a wavelength dispersion in a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, including: means for causing incident light having a known spectral density function S(ω) to be incident on the optical waveguide; means for observing, by the incident light from the means for causing the incident light on the optical waveguide, a signal being proportional to a scattered light amplitude generated at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and a signal being proportional to a scattered light amplitude generated at a second halfway point being different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way; and means for calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by a correlational function between the signal being proportional to the scattered light amplitude at the first halfway point and the signal being proportional to the scattered light amplitude at the second halfway point, the signals being observed by the means for observing the signals which are proportional to the scattered light amplitudes.

In the wavelength dispersion measuring apparatus, the means for observing the signals being proportional to the scattered light amplitudes desirably observes Rayleigh-scattered lights as the scattered lights at the first halfway point and the second halfway point.

In the wavelength dispersion measuring apparatus, the means for observing the signals being proportional to the scattered light amplitudes desirably guides scattered light generated by incident light being incident on the optical waveguide and incident light passing through delay means having delay time equal to propagation time corresponding to twice a length between the incident end of the optical waveguide and the first halfway point to first interference signal detecting means and second interference signal detecting means which have a relative delay difference equal to a distance between the first halfway point and the second halfway point to detect an interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and an interference signal I₂ being proportional to a scattered light amplitude from the second halfway point, and the means for calculating the wavelength dispersion desirably calculates a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by a correlational function between the interference signal I₁ and the interference signal I₂ detected by the means for observing the signals being proportional to the scattered light amplitudes.

In the wavelength dispersion measuring apparatus, variable delay means set on a path of the scattered light or the incident light is arranged, the means for observing the signals being proportional to the scattered light amplitudes desirably respectively detects the plurality of interference signals I₁ and I₂ while causing the variable delay means to change quantity of delay, and the means for calculating the wavelength dispersion desirably calculates a correlational function between the interference signals I₁ and I₂ by calculating the following numerical expression (45) expressing an ensemble average of products of the interference signals I₁ and I₂ detected by the means for observing the signals being proportional to the scattered light amplitudes: I₁I₂*  [Numerical Expression 45]

In the wavelength dispersion measuring apparatus, variable delay means set on a path of the scattered light or the incident light is arranged, the means for observing the signals being proportional to the scattered light amplitudes desirably detects the interference signals I₁ and I₂ as interference signals I while continuously change delay time τ of the variable delay means, and the means for calculating the wavelength dispersion desirably records the interference signals I as a function I(τ) of τ and calculates a correlational function between the interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and the interference signal I₂ being proportional to a scattered light amplitude from the second halfway point by calculating the following numerical expression (47) expressing an ensemble average of the interference signals I₁ and I₂: I₁I₂*  [Numerical Expression 47] by the following numerical expression (46): I ₁ I ₂*=

I(τ_(i))I(τ_(i) +d/c)

_(i)  [Numerical Expression 46] (where < >_(i) expresses an averaging operation about τi) by using the function I(τ).

In the wavelength dispersion measuring apparatus, by using the following numerical expression (48) calculated as the correlational function of the interference signals I₁ and I₂: I₁I₂*  [Numerical Expression 48] a dispersion D in the section extending from the incident end of the optical waveguide to the arbitrary point along the way is desirably calculated by the following numerical expression (49):

$\begin{matrix} {\mspace{79mu}\begin{matrix} {\frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{1}}^{2}}} = \frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{2}}^{2}}}} \\ {= \frac{\int_{- \infty}^{\infty}{{\Gamma(Z)}{\Gamma^{*}\left( {Z - d} \right)}{\mathbb{d}Z}}}{\int_{- \infty}^{\infty}{{{\Gamma(Z)}}^{2}{\mathbb{d}Z}}}} \end{matrix}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 49} \right\rbrack \\ {\mspace{79mu}{{{where},{{\Gamma(z)} \cong {{\exp\left\lbrack {{\mathbb{i}}\;\frac{\omega_{0}}{v_{p}}z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}}}\mspace{79mu}{D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}}}} & \; \end{matrix}$

-   S(ω): spectral density function of known incident light -   d: relative delay difference between first and second interference     signal detecting means -   ω₀: central frequency of incident light -   ν_(g): group velocity of optical waveguide (can be calculated from     refractive index or the like of optical waveguide) -   z₀: length of section extending from incident end of the optical     waveguide to arbitrary point on the way -   λ: central wavelength of incident light     |I₁|², |I₂|²: squares of absolute values of I₁ and I₂

The present invention is a wavelength dispersion measuring method for an optical waveguide capable of measuring a wavelength dispersion in a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, including: a step of causing incident light having a known spectral density function S(ω) to be incident on the optical waveguide; a step of observing, by the incident light being incident on the optical waveguide in the step of causing the incident light to be incident on the optical waveguide, a signal being proportional to a scattered light amplitude generated at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and a signal being proportional to a scattered light amplitude generated at a second halfway point being different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way; and a step of calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by square of correlational function between the signal being proportional to the scattered light amplitude at the first halfway point and the signal being proportional to the scattered light amplitude at the second halfway point, the signals being observed in the step of observing the signals which are proportional to the scattered light amplitudes.

In the wavelength dispersion measuring method, the step of observing the signals being proportional to the scattered light amplitudes is desirably the step of guiding scattered light generated by incident light being incident on the optical waveguide and incident light passing through delay means having delay time equal to propagation time corresponding to twice a length between the incident end of the optical waveguide and the first halfway point to first interference signal detecting means and second interference signal detecting means which have a relative delay difference equal to a distance between the first halfway point and the second halfway point to detect an interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and an interference signal I₂ being proportional to a scattered light amplitude from the second halfway point, and the step of calculating the wavelength dispersion is desirably the step of calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by a correlational function between squares |I₁|² and |I₂|² of absolute values of the interference signals I₁ and I₂.

In the wavelength dispersion measuring method, variable delay means set on a path of the scattered light or the incident light is further arranged, the step of observing the signals being proportional to the scattered light amplitudes is desirably the step of detecting the interference signals I₁ and I₂ as interference signals I while continuously changing delay time τ of the variable delay means, and the step of calculating the wavelength dispersion is desirably the step of recording the interference signals I as a function I(τ) of τ and calculating square of correlational function between the interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and the interference signal I₂ being proportional to a scattered light amplitude from the second halfway point by using the function I(τ).

In the wavelength dispersion measuring method, the correlational function of the squares of the absolute values of the interference signals I₁ and I₂ are calculated to calculate a dispersion D in the section extending from the incident end of the optical waveguide to the arbitrary point along the way by the following numerical expression (50):

$\begin{matrix} \begin{matrix} {\overset{\_}{\left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right)\left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right)} = \overset{\_}{{{I_{1}}^{2}{I_{2}}^{2}} - {I}^{2}}} \\ {= {K\begin{Bmatrix} {\int_{- \infty}^{\infty}{{{\Gamma(z)}} \cdot}} \\ {{{\Gamma\left( {z - d} \right)}\ }{\mathbb{d}z}} \end{Bmatrix}^{2}}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 50} \right\rbrack \\ {{where},} & \; \\ {{{\Gamma(z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack \begin{matrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{matrix}\  \right\rbrack} \right\}{\mathbb{d}\omega}}}}{D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}}} & \; \end{matrix}$

-   S(ω): spectral density function of known incident light -   d: relative delay difference between first and second interference     signal detecting means -   ω₀: central frequency of incident light -   ν_(g): group velocity of optical waveguide (can be calculated from     refractive index or the like of optical waveguide) -   z₀: length of section extending from incident end of the optical     waveguide to arbitrary point along the way -   λ: central wavelength of incident light -   K: constant -   |I₁|², |I₂|²: squares of absolute values of I₁ and I₂.     |I|² :average of |I₁|² and |I₂|²  [Numerical Expression 51]

The present invention is a wavelength dispersion measuring apparatus for an optical waveguide capable of measuring a wavelength dispersion in a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, including: means for causing incident light having a known spectral density function S(ω) to be incident on the optical waveguide; means for observing, by the incident light from the means for causing the incident light on the optical waveguide, a signal being proportional to a scattered light amplitude generated at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and a signal being proportional to a scattered light amplitude generated at a second halfway point being different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way; and means for calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by square of correlational function between the signal being proportional to the scattered light amplitude at the first halfway point and the signal being proportional to the scattered light amplitude at the second halfway point, the signals being observed by the means for observing the signals which are proportional to the scattered light amplitudes.

In the wavelength dispersion measuring apparatus, the means for observing the signals being proportional to the scattered light amplitudes desirably guides scattered light generated by incident light being incident on the optical waveguide and incident light passing through delay means having delay time equal to propagation time corresponding to twice a length between the incident end of the optical waveguide and the first halfway point to first interference signal detecting means and second interference signal detecting means which have a relative delay difference equal to a distance between the first halfway point and the second halfway point to detect an interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and an interference signal I₂ being proportional to a scattered light amplitude from the second halfway point, and the means for calculating the wavelength dispersion desirably measuring a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by a correlational function between squares |I₁|² and |I₂|² of absolute values of the interference signals I₁ and I₂ detected by the means for observing the signals being proportional to the scattered light amplitudes.

In the wavelength dispersion measuring apparatus, variable delay means set on a path of the scattered light or the incident light is further arranged, the means for observing the signals being proportional to the scattered light amplitudes desirably detects the interference signals I₁ and I₂ as interference signals I while continuously changing delay time τ of the variable delay means, and the means for calculating the wavelength dispersion desirably records the interference signals I as a function I(τ) of τ and calculates square of correlational function between the interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and the interference signal I₂ being proportional to a scattered light amplitude from the second halfway point by using the function I(τ).

In the wavelength dispersion measuring apparatus, the means for calculating the wavelength dispersion desirably calculates the correlational function of the squares of the absolute values of the interference signals I₁ and I₂ to calculate a dispersion D in the section extending from the incident end of the optical waveguide to the arbitrary point along the way by the following numerical expression (52):

$\begin{matrix} \begin{matrix} {\overset{\_}{\left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right)\left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right)} = {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} -}} \\ {\overset{\_}{{I}^{2}}} \\ {{= {K\begin{Bmatrix} {\int_{- \infty}^{\infty}{{{\Gamma(z)}} \cdot}} \\ {{{\Gamma\left( {z - d} \right)}}\ {\mathbb{d}z}} \end{Bmatrix}^{2}}}} \end{matrix} & \left\lbrack \;{{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 52} \right\rbrack \\ {{where},} & \; \\ {{{\Gamma(z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}\ {\mathbb{d}\omega}}}}{D = {{- \frac{2\;\pi\; c}{\lambda^{2}}}\beta^{''}}}} & \; \end{matrix}$

-   S(ω): spectral density function of known incident light -   d: relative delay difference between first and second interference     signal detecting means -   ω₀: central frequency of incident light -   ν_(g): group velocity of optical waveguide (can be calculated from     refractive index or the like of optical waveguide) -   z₀: length of section extending from incident end of the optical     waveguide to arbitrary point along the way -   λ: central wavelength of incident light -   K: constant -   |I₁|², |I₂|²: squares of absolute values of I₁ and I₂.     |I|² :average of |I₁|² and |I₂|²  [Numerical Expression 53]

The present invention is a measuring program used in a wavelength dispersion measuring apparatus for an optical waveguide capable of measuring a wavelength dispersion in a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, wherein an interference signal I₁ being proportional to a scattered light amplitude generated by incident light being incident on the optical waveguide at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and an interference signal I₂ being proportional to a scattered light amplitude generated by the incident light being incident on the optical waveguide at a second halfway point different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way are received as numerical values, and the correlational function of the squares of the absolute values of the interference signals I₁ and I₂ to calculate a wavelength dispersion D in the section extending from the incident end of the optical waveguide to the arbitrary point along the way by the following numerical expression (54):

$\begin{matrix} \begin{matrix} {\overset{\_}{\left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right)\left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right)} = {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} -}} \\ {\overset{\_}{{I}^{2}}} \\ {{= {K\begin{Bmatrix} {\int_{- \infty}^{\infty}{{{\Gamma(z)}} \cdot}} \\ {{{\Gamma\left( {z - d} \right)}}\ {\mathbb{d}z}} \end{Bmatrix}^{2}}}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 54} \right\rbrack \\ {{where},} & \; \\ {{{\Gamma(z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} - {\frac{\omega - \omega_{0}}{v_{g}}z}} \right\rbrack} \right\}\ {\mathbb{d}\omega}}}}{D = {{- \frac{2\;\pi\; c}{\lambda^{2}}}\beta^{''}}}} & \; \end{matrix}$

-   S(ω): spectral density function of known incident light -   d: relative delay difference between first and second interference     signal detecting means -   ω₀: central frequency of incident light -   ν_(g): group velocity of optical waveguide (can be calculated from     refractive index or the like of optical waveguide) -   z₀: length of section extending from incident end of the optical     waveguide to arbitrary point along the way -   λ: central wavelength of incident light -   K: constant -   |I₁|², |I₂|²: squares of absolute values of I₁ and I₂.     |I|I² :average of |I₁|² and |I₂|²  [Numerical Expression 55]

The present invention is a measuring program used in a wavelength dispersion measuring apparatus for an optical waveguide capable of measuring a wavelength dispersion in a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, wherein an interference signal I₁ being proportional to a scattered light amplitude generated by incident light being incident on the optical waveguide at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and an interference signal I₂ being proportional to a scattered light amplitude generated by the incident light being incident on the optical waveguide at a second halfway point different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way are received as numerical values, and the correlational function of the interference signals I₁ and I₂ is calculated to calculate a wavelength dispersion D in the section extending from the incident end of the optical waveguide to the arbitrary point along the way by the following numerical expression (56):

$\begin{matrix} \begin{matrix} {\frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{1}}^{2}}} = \frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{2}}^{2}}}} \\ {= \frac{\int_{- \infty}^{\infty}{{\Gamma(Z)}{\Gamma^{*}\left( {Z - d} \right)}{\mathbb{d}Z}}}{\int_{- \infty}^{\infty}{{{\Gamma(Z)}}^{2}{\mathbb{d}Z}}}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 56} \right\rbrack \\ {{where},} & \; \\ {{{\Gamma(z)} \cong {{\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}}z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} - {\frac{\omega - \omega_{0}}{v_{g}}z}} \right\rbrack} \right\}\ {\mathbb{d}\omega}}}}}{D = {{- \frac{2\;\pi\; c}{\lambda^{2}}}\beta^{''}}}} & \; \end{matrix}$

-   S(ω): spectral density function of known incident light -   d: relative delay difference between first and second interference     signal detecting means -   ω₀: central frequency of incident light -   ν_(g): group velocity of optical waveguide (can be calculated from     refractive index or the like of optical waveguide) -   z₀: length of section extending from incident end of the optical     waveguide to arbitrary point along the way -   λ: central wavelength of incident light -   |I₁|², |I₂|²: squares of absolute values of I₁ and I₂.

The present invention is a method capable of measuring a wavelength dispersion of a section extending from an incident end of an optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, including: a step of causing incident light having a known spectral density function S(ω) to be incident on the optical waveguide; a step of observing, by the incident light being incident on the optical waveguide in the step of causing the incident light to be incident on the optical waveguide, a signal being proportional to a scattered light amplitude generated at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and a signal being proportional to a scattered light amplitude generated at a second halfway point being different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way; and a step of calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by square of correlational function between the signal being proportional to the scattered light amplitude at the first halfway point and the signal being proportional to the scattered light amplitude at the second halfway point, the signals being observed in the step of observing the signals which are proportional to the scattered light amplitudes, wherein the step of observing the signals being proportional to the scattered light amplitudes is desirably the step of, when scattered light generated by incident light being incident on the optical waveguide and incident light passing through delay means having delay time equal to propagation time corresponding to twice a length between the incident end of the optical waveguide and the first halfway point are guided to first interference signal detecting means and second interference signal detecting means which have a relative delay difference equal to a distance between the first halfway point and the second halfway point to detect an interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and an interference signal I₂ being proportional to a scattered light amplitude from the second halfway point, detecting the interference signals I₁ and I₂ through a dispersion medium which is inserted on a path on which the scattered lights from the first halfway point and the second halfway point are guided to the first interference signal detecting means and the second interference signal detecting means or a path on which the incident light is guided to the first interference signal detecting means and the second interference signal detecting means in advance and in which a value of the following numerical expression (57): β_(ref)″L_(ref)  [Numerical Expression 57] expressing an accumulated wavelength dispersion is known, and the step of calculating the wavelength dispersion is desirably the step of calculating a correlational function between squares of absolute values of the interference signals I₁ and I₂ to calculate a dispersion D in the section extending from the incident end of the optical waveguides to the arbitrary point along the way by the following Numerical Expression (58):

$\begin{matrix} \begin{matrix} {\overset{\_}{\left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right)\left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right)} = {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} -}} \\ {\overset{\_}{{I}^{2}}} \\ {{= {K\begin{Bmatrix} {\int_{- \infty}^{\infty}{{{\Gamma(z)}} \cdot}} \\ {{{\Gamma\left( {z - d} \right)}}\ {\mathbb{d}z}} \end{Bmatrix}^{2}}}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 58} \right\rbrack \\ {{where},} & \; \\ {{{\Gamma(z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{{\beta^{''}z_{0}} - {\beta_{ref}^{''}L_{ref}}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}\ {\mathbb{d}\omega}}}}{D = {{- \frac{2\;\pi\; c}{\lambda^{2}}}\beta^{''}}}} & \; \end{matrix}$

-   S(ω): spectral density function of known incident light -   d: relative delay difference between first and second interference     signal detecting means -   ω₀: central frequency of incident light -   ν_(g): group velocity of optical waveguide (can be calculated from     refractive index or the like of optical waveguide) -   z₀: length of section extending from incident end of the optical     waveguide to arbitrary point along the way -   λ: central wavelength of incident light -   K: constant -   |I₁|², |I₂|²: squares of absolute values of I₁ and I₂.     |I|² :average of |I₁|² and |I₂|²  [Numerical Expression 59]

The present invention is a wavelength dispersion measuring apparatus for an optical waveguide capable of measuring a wavelength dispersion of a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, including: means for causing incident light having a known spectral density function S(ω) to be incident on the optical waveguide; means for observing, by the incident light from the means for causing the incident light to be incident on the optical waveguide, a signal being proportional to a scattered light amplitude generated at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and a signal being proportional to a scattered light amplitude generated at a second halfway point being different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way; and means for calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by square of correlational function between the signal being proportional to the scattered light amplitude at the first halfway point and the signal being proportional to the scattered light amplitude at the second halfway point, the signals being observed in the means of observing the signals which are proportional to the scattered light amplitudes, wherein the means for observing the signals being proportional to the scattered light amplitudes desirably detects, when scattered light generated by incident light being incident on the optical waveguide and incident light passing through delay means having delay time equal to propagation time corresponding to twice a length between the incident end of the optical waveguide and the first halfway point are guided to first interference signal detecting means and second interference signal detecting means which have a relative delay difference equal to a distance between the first halfway point and the second halfway point to detect an interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and an interference signal I₂ being proportional to a scattered light amplitude from the second halfway point, the interference signals I₁ and I₂ through a dispersion medium which is inserted on a path on which the scattered lights from the first halfway point and the second halfway point are guided to the first interference signal detecting means and the second interference signal detecting means or a path on which the incident light is guided to the first interference signal detecting means and the second interference signal detecting means in advance and in which a value of the following Numerical Expression (60): β_(ref)″L_(ref)  [Numerical Expression 60] expressing an accumulated wavelength dispersion is known, and the means for calculating the wavelength dispersion desirably calculates a correlational function between squares of absolute values of the interference signals I₁ and I₂ to calculate a dispersion D in the section extending from the incident end of the optical waveguides to the arbitrary point along the way by the following Numerical Expression (61):

$\begin{matrix} \begin{matrix} {\overset{\_}{\left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right)\left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right)} = {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} -}} \\ {\overset{\_}{{I}^{2}}} \\ {{= {K\begin{Bmatrix} {\int_{- \infty}^{\infty}{{{\Gamma(z)}} \cdot}} \\ {{{\Gamma\left( {z - d} \right)}}\ {\mathbb{d}z}} \end{Bmatrix}^{2}}}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 61} \right\rbrack \\ {{where},} & \; \\ {{{\Gamma(z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{{\beta^{''}z_{0}} - {\beta_{ref}^{''}L_{ref}}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}\ {\mathbb{d}\omega}}}}{D = {{- \frac{2\;\pi\; c}{\lambda^{2}}}\beta^{''}}}} & \; \end{matrix}$

-   S(ω): spectral density function of known incident light -   d: relative delay difference between first and second interference     signal detecting means -   ω₀: central frequency of incident light -   ν_(g): group velocity of optical waveguide (can be calculated from     refractive index or the like of optical waveguide) -   z₀: length of section extending from incident end of the optical     waveguide to arbitrary point along the way -   λ: central wavelength of incident light -   K: constant -   |I₁|², |I₂|²: squares of absolute values of I₁ and I₂.     |I|² :average of |I₁|² and |I₂|²  [Numerical Expression 62]

A typical example of a scattering phenomenon in an optical waveguide which can be used in the present invention is Rayleigh scattering generated by fluctuation in refractive index in a small space region of the optical waveguide. Furthermore, a scattering phenomenon in which a local scattering coefficient is independent of time, such as a scattering phenomenon generated by incompleteness of, for example, a waveguide structure or a scattering phenomenon generated by containing an impurity can be used in the present invention. In contradiction to this, Brillouin scattering, Raman scattering, or the like known as a scattering phenomenon in an optical waveguide is generated by exciting acoustic phonons or optical phonons in an optical fiber, respectively. However, these local scattering coefficients fluctuate with time because the phonons are repeatedly generated and extinguished in a predetermined life cycle. These phenomena cannot be used in the present invention. This is because a correlational function of a scattered light amplitude or a scattered light intensity is eliminated by a fluctuation in local scattering intensity with time and cannot be measured.

The scattering phenomenon generated by incompleteness of the waveguide structure or the scattering phenomenon generated by containing an impurity decreases with improvement of quality of the waveguide, and cannot be observed in an ideal waveguide. In a high-quality waveguide, these scattering phenomena can be observed at only a limited position of the waveguide. On the other hand, since Rayleigh scattering is always observed at an arbitrary position in an amorphous material, it is assumed that a Rayleigh scattering phenomenon in an optical fiber is used. However, even though another scattering phenomenon is used, these explanations do not essentially change.

In Non-patent Document 2, a technique that analyzes a waveguide having a waveguide dispersion at a high spatial resolution by a configuration similar to that of the present invention is disclosed. In application of the theory disclosed in Non-patent Document 2, when there is a scattering point the position of which can be spatially identified and which is specially stronger than that therearound, a wavelength dispersion of a section till the scattering point can be calculated. On the other hand, in the present invention, by using a scattering phenomenon such as Rayleigh scattering in which respective scattering positions cannot be specified and are uniformly statistically present in a waveguide, a wavelength dispersion in an arbitrary section can be measured by a theory different from that of Non-patent Document 2.

In one of embodiments of the present invention, a variable delay means interlocked with a data acquisition control means is arranged on a path of a the scattered light or wide-band spectral light. The data acquisition control means records outputs from a first interference signal detecting means and a second interference signal detecting means on a data acquisition storing means while changing quantity of delay of the variable delay means.

As the optical waveguide, a slab type or channel type optical waveguide, optical fiber, or the like which is used in a planar lightwave circuit and which confines light in a predetermined medium to propagate the light is conceived.

EFFECT OF THE INVENTION

In a method, an apparatus, and a program for measuring wavelength dispersion of an optical waveguide according to the present invention, a wavelength dispersion in an arbitrary section of an optical waveguide can be measured by a method having an essential principle different from a conventional method. More specifically, there is provided a new method of non-destructively measuring a distribution of wavelength dispersions of an optical waveguide in a direction of length of the optical waveguide. A method, an apparatus, and a program for measuring wavelength dispersion without requiring a high-output light source unlike in a conventional technique can be provided.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is an explanatory block diagram showing a wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 1-1 of the present invention.

FIG. 2 is an explanatory block diagram showing an example of an interference signal detecting means according to Embodiment 1-1 of the present invention.

FIG. 3 is an explanatory block diagram showing a wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 1-2 of the present invention.

FIG. 4 is an explanatory block diagram showing a wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 2-1 of the present invention.

FIG. 5 is an explanatory block diagram showing a wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 2-2 of the present invention.

FIG. 6 is an explanatory block diagram showing a wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 3 of the present invention.

FIG. 7 is a schematic diagram showing a frequency modulated waveform in an OFDR method.

FIG. 8 is a schematic diagram showing a spectral density function in the OFDR method.

REFERENCE NUMERALS

Reference numerals used in the drawings are as follows. 11: measured optical fiber, 12: incident end, 13: halfway point, 14: wide-band spectral light generating means, 15: variable bandpass filter means, 16: optical splitting means, 17: optical circulate means, 18: optical splitting means, 19: first interference signal detecting means, 20: delay means, 21: optical splitting means, 22: variable delay means, 23: second interference signal detecting means, 24: data acquiring/storing means, 25: data calculating means, 26: data acquisition control means, 27: variable delay means, 28, 29: optical splitting/coupling means, 30: interference signal detecting means, 31, 32, 34, 35: photodiodes, 33, 36: balance type photodetectors, 37: movable mirror, 38: stepwise variable delay means, 39: optical phase modulating means, 40: polarization control means, 41: optical amplifying means, 51: absolute-value-square calculating means, 61: narrow-line-width frequency modulated light generating means, 71: reference dispersion giving means, 90 a, b, c, d: reference light, 91 a, b: Rayleigh-scattered light, 101: every-distance Rayleigh-scattered light amplitude analyzing means, 102: every-distance Rayleigh-scattered light amplitude data storing means

BEST MODE FOR CARRYING OUT THE INVENTION

(Embodiment 1-1) Embodiments of the present invention will be described below in detail with reference to the drawings.

FIG. 1 is an explanatory block diagram showing a wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 1-1 of the present invention. This embodiment describes a case in which an optical waveguide is a measured optical fiber 11. However, the optical waveguide is not limited to the measured optical fiber 11, and a slab type or channel type optical waveguide which is used in a planar lightwave circuit and confirms light in a predetermined medium to propagate the light can also be applied.

As shown in FIG. 1, a wavelength dispersion of a section (to be referred to as a measuring section in the embodiment) extending from an incident end 12 of the measured optical fiber 11 to an arbitrary halfway point 13 is measured. Wide-band spectral light (reference light 90 a) having a known spectral density function S(ω) from a wide-band spectral light generating means 14 is incident on the measured optical fiber 11 through a variable bandpass filter means 15, an optical splitting means 16, and an optical circulate means 17.

Rayleigh-scattered light 91 a generated at a halfway point of the measured optical fiber 11 by the wide-band spectral light (reference light 90 b) being incident on the measured optical fiber 11 is input to a first interference signal detecting means 19 through the optical circulate means 17 and an optical splitting means 18. Wide-band spectral light (reference light 90 c) split by the optical splitting means 16 is input to the first interference signal detecting means 19 through a delay means 20 having delay time equal to propagation time corresponding to twice a length (measuring section length) from the incident end 12 of the measured optical fiber 11, a variable delay means 27 interlocked with a data acquisition control means 26, and an optical splitting means 21. The Rayleigh-scattered light 91 b split by the optical splitting means 18 is input to a second interference signal detecting means 23 through a variable delay means 22. A wide-band spectral light (reference light 90 d) split by the optical splitting means 21 is input to the second interference signal detecting means 23.

In the first interference signal detecting means 19 and the second interference signal detecting means 23, signals being proportional to Rayleigh-scattered light amplitudes are detected to be input to a data acquiring/storing means 24. In the data acquiring/storing means 24 and a data calculating means 25, a wavelength dispersion of the section extending from the incident end 12 of the measured optical fiber 11 to the arbitrary halfway point 13 is calculated by a correlational function of signals being proportional to the Rayleigh-scattered light amplitudes generated at a plurality of halfway points. The variable delay means 27 and the data acquiring/storing means 24 are controlled by the data acquisition control means 26. The data acquisition control means 26 records outputs from the first interference signal detecting means 19 and the second interference signal detecting means 23 on the data acquiring/storing means 24 while changing quantity of delay of the variable delay means 27.

The first interference signal and the second interference signal detected by the first interference signal detecting means 19 and the second interference signal detecting means 23 are stored by the data acquiring/storing means 24 and guided to the data calculating means 25. In the data calculating means 25, a correlational function between the first interference signal and the second interference signal is calculated by a numerical operation to make it possible to calculate a wavelength dispersion of the measured optical fiber 11 of a predetermined section on the basis of a theoretical expression (will be described later). When the first interference signal and the second interference signal are represented by I₁ and I₂, respectively, the correlational function between the first interference signal and the second interference signal is given by the following Numerical Expression (80) expressing an ensemble average of products thereof: I₁I₂*  [Numerical Expression 80] where ⋆ is a complex conjugate sign. In order to calculate the ensemble average, a large number of interference signals must be acquired as data while changing quantity of delay as described above.

According to the theory which is disclosed in Embodiment 1-1 for the first time, the correlational function between the first interference signal I₁ and the second interference signal I₂ is associated with a dispersion D in the measuring section by the following relational expression (Numerical Expression (81)):

$\begin{matrix} {\begin{matrix} {\frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{1}}^{2}}} = \frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{2}}^{2}}}} \\ {= \frac{\int_{- \infty}^{\infty}{{\Gamma(Z)}{\Gamma^{*}\left( {Z - d} \right)}{\mathbb{d}Z}}}{\int_{- \infty}^{\infty}{{{\Gamma(Z)}}^{2}{\mathbb{d}Z}}}} \end{matrix}{{where},{{\Gamma(z)} \cong {{\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}}z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} - {\frac{\omega - \omega_{0}}{v_{g}}z}} \right\rbrack} \right\}\ {\mathbb{d}\omega}}}}}}{D = {{- \frac{2\;\pi\; c}{\lambda^{2}}}\beta^{''}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 81} \right\rbrack \end{matrix}$

-   S(ω): spectral density function of known wide-band spectral light -   d: relative delay difference between first and second interference     signal detecting means 19 and 23 -   ω₀: central frequency of wide-band spectral light -   ν_(g): group velocity of the measured optical fiber 11 (can be     calculated from refractive index or the like of the measured optical     fiber 11) -   z₀: length of measuring section extending from the incident end 12     of the measured optical fiber 11 -   λ: central wavelength of wide-band spectral light -   |I₁|², |I₂|²: squares of absolute values of I₁ and I₂

The wide-band spectral light generating means 14 is realized by, for example, a super-luminescent diode and generates a wide-band spectral continuous light having a spectral width of about several ten nanometers. The variable bandpass filter means 15 cuts out spectral light having a central frequency ω₀ and a bandwidth Δ_(ω). The central frequency ω₀ of the spectral light can be varied by the data acquisition control means 26. It is assumed that the spectral density function S(ω) of the continuous light cut out by the variable bandpass filter means 15 is known in advance. In general, an output from the variable bandpass filter means 15 can be approximated by a Gaussian profile and represented as the following Numerical Expression (82).

$\begin{matrix} {{S(\omega)} = \frac{\exp\left\lbrack \frac{- \left( {\omega - \omega_{0}} \right)^{2}}{2\left( {\Delta\;\omega} \right)^{2}} \right\rbrack}{\sqrt{2\;\pi}\Delta\;\omega}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 82} \right\rbrack \end{matrix}$

Here, in Numerical Expression (82), reference symbol ω₀ is a central frequency which can be controlled by the data acquisition control means 26. Reference symbol Δ_(ω) is a bandwidth.

A signal being proportional to a Rayleigh-scattered light amplitude generated at a first halfway point in the measuring section of the measured optical fiber 11 by the wide-band spectral light (reference light 90 b) being incident on the measured optical fiber 11 and a signal being proportional to a Rayleigh-scattered light amplitude generated at a second halfway point included in the measuring section of the measured optical fiber 11 and being different from the first halfway point are observed. In order to realize this, in Embodiment 1-1, the Rayleigh-scattered light 91 a generated by the wide-band spectral light (reference light 90 b) being incident on the measured optical fiber 11 and the wide-band spectral light (reference light 90 c) transmitted through the delay means 20 having delay time equal to propagation time corresponding to twice the measuring section and the variable delay means 27 are guided to the first interference signal detecting means 19 and the second interference signal detecting means 23 which have a relative optical path length difference being equal to a distance between the first halfway point and the second halfway point to detect signals as the interference signal I₁ being proportional to the Rayleigh-scattered light amplitude from the first halfway point and the interference signal I₂ being proportional to the Rayleigh-scattered light amplitude from the second halfway point.

More specifically, after the wide-band spectral light (reference light 90 a) passes through the variable bandpass filter means 15, the wide-band spectral light (reference light 90 a) is split by two by means of the optical splitting means 16. One of the split lights is guided to the measured optical fiber 11. At this time, the Rayleigh-scattered light 91 a caused by a fluctuation in refractive index of glass at an arbitrary point of the measured optical fiber 11 is generated, and the generated Rayleigh-scattered light 91 a is guided to the first and second interference signal detecting means 19 and 23 through the optical circulate means 17. The other wide-band spectral light (reference light 90 c) of the lights split by the optical splitting means 16 is guided to the first and second interference signal detecting means 19 and 23 through the delay means 20 having delay time equal to propagation time corresponding to twice the measuring section length and the variable delay means 27. The delay means 20 having delay time equal to the propagation time corresponding to twice the measuring section length is set so as to make the optical path lengths of the following optical paths 1 and 2 equal to each other.

Optical path 1: optical splitting means 16→optical circulate means 17→arbitrary reflection point (of Rayleigh scattering) of measured optical fiber 11→optical circulate means 17→first and second interference signal detecting means 19 and 23. Optical path 2: optical splitting means 16→first and second interference signal detecting means 19 and 23.

More specifically, this can be realized by changing connections of a fixed delay line of an optical fiber. The length of the optical fiber is changed to make it possible to arbitrarily change a measuring section.

The reference light 90 c and the Rayleigh-scattered light 91 a are split by two by the optical splitting means 21 and 18, respectively. The first interference signal detecting means 19 and the second interference signal detecting means 23 which have a relative optical path length difference equal to a distance between the first halfway point and the second halfway point of the measured optical fiber 11 detect interference signals generated by the reference lights 90 c and 90 d and the Rayleigh-scattered lights 91 a and 91 b. When the delay means 20 having delay time equal to propagation time corresponding to twice the measuring section length and the variable delay means 27 are inserted into the path of the reference light 90 c, the interference signal I₁ being proportional to the Rayleigh-scattered light amplitude from the first halfway point of the measured optical fiber 11 is output from the first interference signal detecting means 19. On the other hand, from the second interference signal detecting means 23, the interference signal I₂ being proportional to the Rayleigh-scattered light amplitude is output from the second halfway point of the measured optical fiber 11.

In this embodiment, the interference signals detected by the first interference signal detecting means 19 and the second interference signal detecting means 23 are finally guided to the data calculating means 25 to make it possible to derive a wavelength dispersion of a predetermined section of the measured optical fiber 11 by a theoretical expression given by a knowledge (will be described later).

The first and second interference signal detecting means 19 and 23 can be constituted by a heterodyne detecting technique or a homodyne detecting technique which is a known technique. As constituent elements thereof, several variations can be conceived.

FIG. 2 is an explanatory block diagram showing an example of an interference signal detecting means according to Embodiment 1-1 of the present invention. The same reference numerals as in FIG. 1 denote the same parts in FIG. 2, and a description thereof will not be repeated. In FIG. 2, as the first interference signal detecting means 19 (FIG. 1), a balance type photodetector 33 constituted by photodiodes 31 and 32 is arranged. As the second interference signal detecting means 23 (in FIG. 1), a balance type photodetector 36 constituted by photodiodes 34 and 35 is arranged. As a variable delay means 22 (FIG. 1), a movable mirror (delay giving unit) 37 is arranged. As the delay means 20 (FIG. 1) having delay time equal to propagation time corresponding to twice the length of the measuring section extending from the incident end 12 of the measured optical fiber 11, a stepwise variable means 38 is arranged. Between the stepwise variable delay means 38 and the optical splitting means 21, the variable delay means 27, an optical phase modulating means (ω_(m)) 39, and a polarization control means 40 are arranged. An optical amplifying means 41 is arranged between the wide-band spectral light generating means 14 and the variable bandpass filter means 15.

As shown in FIG. 2, the optical phase modulating means (ω_(m)) 39 phase-modulates the reference light 90 c by a sine wave consisting of a specific frequency ω_(m). The phase-modulated lights are coupled by optical splitting/coupling means 28 and 29 and photoelectrically converted by the balance type photodetector 33 and a balance type photodetector 36. Thereafter, the intensity of a frequency component consisting of the frequency ω_(m) in the data acquiring/storing means 24 and the data calculating means 25 is measured by a spectral analyzing means (not shown) to make it possible to observe coherent correlation. The polarization control means 40 is to cope with a fluctuation in polarization caused by an optical fiber and adjusts a polarization state of the reference light 90 c to such that the interference signal intensity is maximum.

The distance between the first halfway point and the second halfway point of the measured optical fiber 11, i.e., a relative delay difference d between the first interference signal detecting means 19 and the second interference signal detecting means 23 is preferably set as follows in accordance with a bandwidth Δ_(ω) of a wide-band spectral light to be used by changing delay time τ in the variable delay means 22. Even though these parameters are set at any values, measurement is not impossible. However, there is a setting range in which the measurement can be performed at relatively high accuracy. A measuring section length is represented by z₀, and a predicted wavelength dispersion in the measuring section is given by the following Numerical Expression (83). {tilde over (D)}  [Numerical Expression 83] When a spectral density function is assumed as Numerical Expression (82), observation with the maximum accuracy can be performed when the following Numerical Expression (84) representing coherence time of wide-band spectral light when the wide-band spectral light reciprocates in a measuring section is equal to d/c (c is velocity of light). This is because observation at high accuracy can be performed since quantity of observation given by Numerical Expression (89) (will be described later) is 1 when d=0, is 0 when d=∞, and is the maximum differential coefficient to the wavelength dispersion when d/c=T_(c).

$\begin{matrix} {T_{c} = {\frac{1}{\Delta\;\omega}\left\lbrack {1 + \left( \frac{2\;\ln\; 2\;\lambda^{2}{\Delta\omega}^{2}\overset{\sim}{D}z_{0}}{\pi\; c} \right)^{2}} \right\rbrack}^{\frac{1}{2}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 84} \right\rbrack \end{matrix}$

-   λ: central wavelength of incident light -   z₀: length of measuring section extending from the incident end 12     of the measured optical fiber 11

A set range required in Numerical Expression (84) has a common width. When the set range is approximately given by 0.1 T_(c)<d/c<10 T_(c), measurement can be sufficiently performed. In actual measurement, when the wavelength dispersion cannot be predicted at all, measurement may be performed any number of times while changing d or Δ_(ω). When, as numerical examples, Δ_(ω)=100 GHz/2π, λ=1.55 μm, and L=1 km are given, and when the Numerical Expression (83) is given by a condition expressed by the following expression (85), d is about 25 mm. {tilde over (D)}=1 ps/nm/km  [Numerical Expression 85]

In Embodiment 1-1, a wavelength dispersion of a section extending from the incident end 12 of the optical fiber to the arbitrary halfway point 13 is calculated from the following Numerical Expression (86) expressing a correlational function between the interference signal I₁ being proportional to the Rayleigh-scattered light amplitude generated at the first halfway point and the interference signal I₂ being proportional to the Rayleigh-scattered light amplitude generated at the second halfway point. I₁I₂*  [Numerical Expression 86]

In order to calculate the correlational function between the interference signal I₁ being proportional to the Rayleigh-scattered light amplitude generated at the first halfway point and the interference signal I₂ being proportional to the Rayleigh-scattered light amplitude generated at the second halfway point, the variable delay means 27 is arranged on the path of the scattered light or the wide-band spectral light, the plurality of interference signals I₁ and I₂ are detected while changing quantity of delay of the variable delay means 27, and an ensemble averages of these signals is calculated. In this manner, the correlational function between the interference signals I₁ and I₂ can be calculated as the following Numerical Expression (87). I₁I₂*  [Numerical Expression 87]

For this reason, the variable delay means 27 can be controlled by the data acquisition control means 26. The quantity of delay of the variable delay means 27 for each measurement must be larger than a coherent length of the wide-band spectral light because independent Rayleigh scattering must be observed every measurement. This delay distance is about 18 mm (corresponding value in vacuum) when, for example, the spectral width of the wide-band spectral light is Δ_(ω)=100 GHz/2π.

The interference signals I₁ and I₂ observed in Embodiment 1-1 are complex variables having amplitudes and phases. According to the configuration in FIG. 2, although a current observed by the balance type photodetectors 33 and 36 is vibrated at a frequency ω_(m), the amplitude (intensity) and the phase of the current are observed to determine the interference signals I₁ and I₂ which are complex variables.

In Embodiment 1-1, from Numerical Expression (88) expressing the correlational function calculated as described above, the dispersion D in the measuring section can be calculated by the following theoretical expression Numerical Expression (89):

$\begin{matrix} \overset{\_}{I_{1}I_{2}^{*}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 88} \right\rbrack \\ {{\frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{1}}^{2}}} = {\frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{2}}^{2}}} = \frac{\int_{- \infty}^{\infty}{{\Gamma(Z)}{\Gamma^{*}\left( {Z - d} \right)}{\mathbb{d}Z}}}{\int_{- \infty}^{\infty}\left| {\Gamma(Z)} \middle| {}_{2}{\mathbb{d}Z} \right.}}}{{where},{{\Gamma(z)} \cong {{\exp\left\lbrack {{\mathbb{i}}\;\frac{\omega_{0}}{v_{p}}z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}}}{D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 89} \right\rbrack \end{matrix}$

-   S(ω): spectral density function of known wide-band spectral light -   d: relative delay difference between first and second interference     signal detecting means 19 and 23 -   ω₀: central frequency of wide-band spectral light -   ν_(g): group velocity of the measured optical fiber 11 (can be     calculated from refractive index or the like of the measured optical     fiber 11) -   z₀: length of measuring section extending from the incident end 12     of the measured optical fiber 11 -   λ: central wavelength of wide-band spectral light -   |I₁|², I₂|²: squares of absolute values of I₁ and I₂

The reason why the above theoretical expression is appropriate will be described below. The waveform of the wide-band spectral light can be represented as the following Numerical Expression (90). Ψ(t)=∫_(−∞) ^(∞) A(ω)exp{i[k(ω)z−ωt]}dω  [Numerical Expression 90]

In this expression, reference symbol A(ω) denotes an electric field amplitude component which is Fourier-dissolved. The following Numerical Expression (91) is established.

$\begin{matrix} {{S(\omega)} = \frac{{{A(\omega)}}^{2}}{2\pi}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 91} \right\rbrack \end{matrix}$

In Numerical Expression (90), reference symbol k is a propagation constant, and the path of the reference light 90 c is nondispersive. For this reason, the following Numerical Expression (92) is established.

$\begin{matrix} {{k(\omega)} = {\frac{\omega_{0}}{v_{p}} + \frac{\omega - \omega_{0}}{v_{g}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 92} \right\rbrack \end{matrix}$

The Rayleigh-scattered lights 91 a from respective points z₁ of the optical fiber is expressed by the following Numerical Expression (93) in the inputs to the first and second interference signal detecting means 19 and 23: Ψ_(r) (t)=α(z _(i))·∫_(−∞) ^(∞) A(ω)exp{i[2β(ω)z _(i) −ωt]}dω  [Numerical Expression 93] In this equation, a(z_(i)) denotes a scattering coefficient (local refractive index) in z₁. In Numerical Expression (93), β(ω) denotes a propagation constant of the optical fiber. According to the Taylor expansion, β(ω) can be represented as the following Numerical Expression (94):

$\begin{matrix} {{\beta(\omega)} = {\frac{\omega_{0}}{v_{p}} + \frac{\omega - \omega_{0}}{v_{g}} + {\frac{\beta^{''}}{2}\left( {\omega - \omega_{0}} \right)^{2}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 94} \right\rbrack \end{matrix}$

In Numerical Expression (94), reference symbol v_(p) denotes a phase velocity of light, reference numeral ν_(g) denotes a group velocity of light, and a value β″ has a relationship given by the following Numerical Expression (95) to the wavelength dispersion D. When the value β″ is calculated, the D can be calculated:

$\begin{matrix} {D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 95} \right\rbrack \end{matrix}$

When a total sum of Numerical Expression (93) with respect to z_(i) is calculated, an amplitude of the Rayleigh-scattered light 91 a is obtained as follows:

$\begin{matrix} {{\Psi_{total}(t)} = {\sum\limits_{i}{{a\left( z_{i} \right)} \cdot {\int_{- \infty}^{\infty}{{A(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{2{\beta(\omega)}z_{i}} -} \\ {\omega\; t} \end{bmatrix}} \right\}\ {\mathbb{d}\omega}}}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 96} \right\rbrack \end{matrix}$ Outputs from the first interference signal detecting means and the second interference signal detecting means 23 are given by the following Numerical Expression (97): I ₁=

Ψ_(total)(t)+Ψ(t)|²|

I ₂=

Ψ_(total)(t)+Ψ(t+d/c)|²

  [Numerical Expression 97] (d expresses delay length obtained by the variable delay means 22, c expresses velocity of light, and < > expresses time averaging operation.). For this reason, when Numerical Expression (90) and Numerical Expression (96) are assigned to Numerical Expression (97) to neglect a DC component, the following Numerical Expression 98 is given:

$\begin{matrix} \begin{matrix} \begin{matrix} {I_{1} = \left\langle {\begin{matrix} {{\sum\limits_{i}{{a\left( z_{i} \right)} \cdot {\int_{- \infty}^{\infty}{{A(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{2{\beta(\omega)}z_{i}} -} \\ {\omega\; t} \end{bmatrix}} \right\}\ {\mathbb{d}\omega}}}}} +} \\ {\int_{- \infty}^{\infty}{{A(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack {{{k(\omega)}z_{0}} - {\omega\; t}} \right\rbrack} \right\}\ {\mathbb{d}\omega}}} \end{matrix}}^{2} \right\rangle} \\ {= {\sum\limits_{i}{{a\left( z_{i} \right)}{\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}}{2}\left( {\omega - \omega_{0}} \right)^{2}z_{i}} +} \\ \begin{matrix} {{\frac{\omega - \omega_{0}}{v_{g}}z_{i}} +} \\ {\frac{\omega_{0}}{v_{p}}z_{i}} \end{matrix} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}}} \\ {= {\sum\limits_{i}{{a\left( z_{i} \right)}{\Gamma\left( z_{i} \right)}}}} \\ {I_{2} = \left\langle {\begin{matrix} {{\sum\limits_{i}{{a\left( z_{i} \right)} \cdot {\int_{- \infty}^{\infty}{{A(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack {{2{\beta(\omega)}z_{i}} - {\omega\; t}} \right\rbrack} \right\}\ {\mathbb{d}\omega}}}}} +} \\ {\int_{- \infty}^{\infty}{{A(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack {{{k(\omega)}\left( {z + d} \right)} - {\omega\; t}} \right\rbrack} \right\}\ {\mathbb{d}\omega}}} \end{matrix}}^{2} \right\rangle} \\ {= {\sum\limits_{i}{{a\left( z_{i} \right)}{\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}}{2}\left( {\omega - \omega_{0}} \right)^{2}\left( {z_{i} - d} \right)} +} \\ \begin{matrix} {{\frac{\omega - \omega_{0}}{v_{g}}\left( {z_{i} - d} \right)} +} \\ {\frac{\omega_{0}}{v_{p}}\left( {z_{i} - d} \right)} \end{matrix} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}}} \\ {= {\sum\limits_{i}{{a\left( z_{i} \right)}{\Gamma\left( {z_{i} - d} \right)}{\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}}\left( {z_{i} - d} \right)} \right\rbrack}}}} \end{matrix} \\ {{where},{{\Gamma(z)} = {{\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}}z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}}{2}\left( {\omega - \omega_{0}} \right)^{2}z} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}}} \end{matrix} & \begin{bmatrix} \begin{matrix} {Numerical} \\ {Expression} \end{matrix} \\ 98 \end{bmatrix} \end{matrix}$

In general, Γ(z) is called a coherence function of wide-band spectral light. In this case, when it is assumed that β″z/2 changes by only quantity which is considerably smaller than 2π with reference to a change of z which is almost equal to d, Γ(z) may be expressed by the following numerical expression:

$\begin{matrix} {{\Gamma(z)} \cong {{\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}}z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 99} \right\rbrack \end{matrix}$

A correlational function of these interference signals is given by the following numerical expression:

$\begin{matrix} \begin{matrix} {\overset{\_}{I_{1}I_{2}^{*}} = \overset{\_}{\sum\limits_{i}{{a\left( z_{i} \right)}{\Gamma\left( z_{i} \right)} \times {\sum\limits_{j}{{a^{*}\left( z_{j} \right)}{\Gamma^{*}\left( {z_{j} - d} \right)}}}}}} \\ {= {\overset{\_}{\sum\limits_{i}{{{a\left( z_{i} \right)}}^{2}{\Gamma\left( z_{i} \right)}{\Gamma^{*}\left( {z_{i} - d} \right)}}} +}} \\ {\overset{\_}{\sum\limits_{i \neq j}{{a\left( z_{i} \right)}{a^{*}\left( z_{j} \right)}{\Gamma\left( z_{i} \right)}{\Gamma^{*}\left( {z_{i} - d} \right)}}}} \\ {= {{\sum\limits_{i}{\overset{\_}{a^{2}\left( z_{i} \right)}{\Gamma\left( z_{i} \right)}{\Gamma^{*}\left( {z_{i} - d} \right)}}} +}} \\ {\sum\limits_{i \neq j}{\overset{\_}{{a\left( z_{i} \right)}{a^{*}\left( z_{j} \right)}}{\Gamma\left( z_{i} \right)}{\Gamma^{*}\left( {z_{i} - d} \right)}}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 100} \right\rbrack \end{matrix}$ Since α(z_(i))α*=(z_(j)) α(z_(i))α*=(z_(j))=0 is satisfied due to randomness of a(z_(i)), the following equation is satisfied:

$\overset{\_}{I_{1}I_{2}^{*}} = {\sum\limits_{i}{\overset{\_}{{{a\left( z_{i} \right)}}^{2}}{\Gamma\left( z_{i} \right)}{\Gamma\left( {z_{i} - d} \right)}}}$ When a standard deviation of a(z₁) is represented by σ, the following numerical expression is satisfied: |α(z _(i))|² =σ²  [Numerical Expression 101] When the following numerical expression is replaced with integration,

$\begin{matrix} \sum\limits_{i} & {\left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{11mu} 102} \right\rbrack\;} \end{matrix}$ the following numerical expression is obtained: I ₁ I ₂*=σ²∫_(−∞) ^(∞Γ)(z)Γ*(z−d)dZ  [Numerical Expression 103] When d=0 is given in Numerical Expression (103), the following numerical expression is satisfied:

$\begin{matrix} {\overset{\_}{{I_{1}}^{2}} = {\overset{\_}{{I_{2}}^{2}} = {{\sigma^{2}{\int_{- \infty}^{\infty}{{{\Gamma(z)}}^{2}\ {\mathbb{d}Z}}}} = {\sigma^{2}{\int_{- \infty}^{\infty}{{{\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} \begin{matrix} \frac{\beta^{''}z_{0}}{2} \\ {\left( {\omega - \omega_{0}} \right)^{2} -} \end{matrix} \\ {\frac{\omega - \omega_{0}}{v_{g}}\ z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}^{2}\ {\mathbb{d}Z}}}}}}} & {\left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 104} \right\rbrack\;} \end{matrix}$ According to Numerical Expression (103) and Numerical Expression (104), the following numerical expression is satisfied:

$\quad\begin{matrix} \begin{matrix} {\frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{1}}^{2}}} = \frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{1}}^{2}}}} \\ {= \frac{\int_{- \infty}^{\infty}{{\Gamma(z)}{\Gamma^{*}\left( {z - d} \right)}\ {\mathbb{d}Z}}}{\int_{- \infty}^{\infty}{{{\Gamma(z)}}^{2}\ {\mathbb{d}Z}}}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 105} \right\rbrack \end{matrix}$ When Γ(z) is restored according to Numerical Expression (98) of Γ(z), the following numerical expression is obtained:

$\begin{matrix} {\frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{1}}^{2}}} = \frac{\begin{matrix} {\int_{- \infty}^{\infty}{\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}\ {\mathbb{d}\omega}}}} \\ {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}\ {\mathbb{d}\omega}{\mathbb{d}Z}}} \end{matrix}}{\int_{- \infty}^{\infty}{{\ {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} \begin{matrix} \frac{\beta^{''}z_{0}}{2} \\ {\left( {\omega - \omega_{0}} \right)^{2} -} \end{matrix} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}\ {\mathbb{d}\omega}}}}^{2}{\mathbb{d} Z}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 106} \right\rbrack \end{matrix}$ Numerical Expression (106) is an integral equation to β″ when S(ω), z₀, and dare known. According to a numerical calculation, the equation can be solved with respect to β″. Therefore, according to the means described above, quantity given by the following expression:

$\begin{matrix} \frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{1}}^{2}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 107} \right\rbrack \end{matrix}$ is experimentally calculated. The quantity is assigned to Numerical Expression (107) to solve Numerical Expression (106) with respect to β″, so that a dispersion β″ of a section till z₀ can be calculated.

(Embodiment 1-2) FIG. 3 is an explanatory block diagram showing a wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 1-2 of the present invention. This embodiment will describe a case in which an optical waveguide is a measured optical fiber 11. The optical waveguide is not limited to the measured optical fiber 11, and a slab type or channel type optical waveguide which is used in a planar lightwave circuit and which confines light in a predetermined medium to propagate the light can also be applied. The same reference numerals as in Embodiment 1-1 denote the same parts in Embodiment 1-2.

As shown in FIG. 3, a wavelength dispersion of a section (to be referred as a measuring section hereinafter in the embodiment) extending from an incident end 12 of a measured optical fiber 11 to an arbitrary halfway point 13 is measured. Wide-band spectral light (reference light 90 a) having a known spectral density function S(ω) from a wide-band spectral light generating means 14 is incident on the measured optical fiber 11 through a variable bandpass filter means 15, an optical splitting means 16, and an optical circulate means 17. Rayleigh-scattered light 91 a generated at a halfway point of the measured optical fiber 11 by the wide-band spectral light (reference light 90 b) being incident on the measured optical fiber 11 is input to an interference signal detecting means 30 serving as a first interference signal detecting means and a second interference signal detecting means through the optical circulate means 17. Wide-band spectral light (reference light 90 c) split by the optical splitting means 16 is input to the interference signal detecting means 30 through a delay means 20 having delay time equal to propagation time corresponding to twice a length of a measuring section from the incident end 12 of the measured optical fiber 11 and a variable delay means 27 interlocked with a data acquisition control means 26.

In the interference signal detecting means 30, signals being proportional to Rayleigh-scattered light amplitudes are detected to be input to a data acquiring/storing means 24 and a data calculating means 25. In the data acquiring/storing means 24 and the data calculating means 25, a wavelength dispersion of the section extending from the incident end 12 of the measured optical fiber 11 to the arbitrary halfway point 13 is calculated by a correlational function of signals being proportional to the Rayleigh-scattered light amplitudes generated at a plurality of halfway points. The variable delay means 27 and the data acquiring/storing means 24 are controlled by the data acquisition control means 26. The data acquisition control means 26 records outputs from the interference signal detecting means 30 on the data acquiring/storing means 24 while changing quantity of delay of the variable delay means 27. The interference signals detected by the interference signal detecting means 30 are stored by the data acquiring/storing means 24 and guided to the data calculating means 25. In the data calculating means 25, a correlational function between the interference signals is calculated by a numerical operation to calculate a dispersion of the measured optical fiber 11 in the predetermined section.

The interference signal detecting means 30 detects the interference signal I in each case while the variable delay means 27 interlocked with the data acquisition control means 26 continuously changes delay time τ to cause the data acquiring/storing means 24 to record the interference signals I as functions of delay time τ (=d/c). At this time, the interference signals I are periodically oscillated at intervals of wide-band spectral light (reference light 90 a), and an amplitude and a phase are measured to observe a complex number I. Therefore, in case of Embodiment 1-2, as accuracy of the variable delay means 27 interlocked with the data acquisition control means 26, accuracy higher than that of a wavelength of light is required. However, this accuracy can be sufficiently achieved by a step pulse motor or the like. In this manner, a plurality of numerical values at two arbitrary points time d/c distant from each other are extracted from I(τ) recorded as functions of the delay time τ. When an ensemble average of these numerical values is calculated, a correlational function between the values time d/c distant from each other in the above Embodiment 1-1: I₁I₂*  [Numerical Expression 108] can be calculated. More specifically, by using the measured I(τ), by the following numerical expression: I ₁ I ₂*=

I(τ_(i))I(τ_(i) +d/c)

_(i)  [Numerical Expression 109] (where < >_(i) is an averaging operation with respect to τ_(j)), the following numerical expression: I₁I₂*  [Numerical Expression 110] may be calculated. A theoretical calculation performed after a value given by the following expression: I₁I₂*  [Numerical Expression 111] is the same as that in the above-described Embodiment 1-1.

The present invention is not directly limited to the example of the above embodiment. In an execution phase, the invention can be embedded by modifying constituent elements within range without departing from the spirit and scope of the invention. Appropriate combinations of a plurality of constituent elements disclosed in the above embodiment can constitute various inventions. For example, several constituent elements may be omitted from all the constituent elements disclosed in the embodiments. Furthermore, the constituent elements across different embodiments may be appropriately combined to each other.

(Embodiment 2-1) An embodiment of the present invention will be described below with reference to the drawings.

FIG. 4 is an explanatory block diagram showing a wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 2-1 of the present invention. This embodiment describes a case in which an optical waveguide is a measured optical fiber 11. However, the optical waveguide is not limited to the measured optical fiber 11, and a slab type or channel type optical waveguide which is used in a planar lightwave circuit and confirms light in a predetermined medium to propagate the light can also be applied.

As shown in FIG. 4, a wavelength dispersion of a section (to be referred to as a measuring section in the embodiment) extending from an incident end 12 of the measured optical fiber 11 to an arbitrary halfway point 13 is measured. Wide-band spectral light (reference light 90 a) having a known spectral density function S(ω) from a wide-band spectral light generating means 14 is incident on the measured optical fiber 11 through a variable bandpass filter means 15, an optical splitting means 16, and an optical circulate means 17. Rayleigh-scattered light 91 a generated at a halfway point of the measured optical fiber 11 by the wide-band spectral light (reference light 90 b) being incident on the measured optical fiber 11 is input to a first interference signal detecting means 19 through the optical circulate means 17 and an optical splitting means 18. Wide-band spectral light (reference light 90 c) split by the optical splitting means 16 is input to the first interference signal detecting means 19 through a delay means 20 having delay time equal to propagation time corresponding to twice a length (measuring section length) from the incident end 12 of the measured optical fiber 11, a variable delay means 27 interlocked with a data acquisition control means 26, and an optical splitting means 21. The Rayleigh-scattered light 91 b split by the optical splitting means 18 is input to a second interference signal detecting means 23. A wide-band spectral light (reference light 90 d) split by the optical splitting means 21 is input to the second interference signal detecting means 23 through a variable delay means 22.

In the first interference signal detecting means 19 and the second interference signal detecting means 23, signals being proportional to Rayleigh-scattered light amplitudes are detected to be input to a absolute-value-square calculating means 51, a data acquiring/storing means 24, and a data calculating means 25. In the absolute-value-square calculating means 51, the data acquiring/storing means 24, and the data calculating means 25, a wavelength dispersion of the section extending from the incident end 12 of the measured optical fiber 11 to the arbitrary halfway point 13 is calculated by a correlational function of signals being proportional to the Rayleigh-scattered light intensities generated at a plurality of halfway points. The variable delay means 27 and the data acquiring/storing means 24 are controlled by the data acquisition control means 26. The data acquisition control means 26 records outputs from the first interference signal detecting means 19 on the data acquiring/storing means 24 while changing quantity of delay of the variable delay means 27.

The first interference signal and the second interference signal detected by the first interference signal detecting means 19 and the second interference signal detecting means 23 are stored by the data acquiring/storing means 24. In the absolute-value-square calculating means 51 and the data calculating means 25, a correlational function between the first interference signal and the second interference signal is calculated by a numerical operation to make it possible to calculate a dispersion of the measured optical fiber 11 of a predetermined section on the basis of a theoretical expression (will be described later).

In the wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 2-1 of the present invention it is characterized that a first interference signal |I₂|² and a second interference signal |I₂|² detected at this time are proportional to squares of absolute values of the Rayleigh-scattered light amplitudes. This merely means that squares of absolute values of currents detected by the first and second interference signal detecting means 19 and 23 may be measured by the absolute-value-square calculating means 51. Therefore, in the wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 2-1 of the present invention, the first interference signal |I₁|² and the second interference signal |I₂|² are positive real numbers.

In the wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 2-1 of the present invention, a dispersion D of an arbitrary section is calculated by a correlational function between the first interference signal |I₁|² and the second interference signal |I₂|² obtained by the means described above using the following relational expressions (numerical expressions (112) and (113)).

$\begin{matrix} \left\lbrack {{Numerical}{\mspace{11mu}\;}{Expression}\mspace{14mu} 112} \right\rbrack & \; \\ {\begin{matrix} {\overset{\_}{\left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right)\left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right)} = {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} - \overset{\_}{{I}^{2}}}} \\ {= {K\left\{ {\int_{- \infty}^{\infty}{{{{\Gamma(z)}} \cdot {{\Gamma\left( {z - d} \right)}}}\ {\mathbb{d}z}}} \right\}^{2}}} \end{matrix}{{where},}} & (A) \\ {{\Gamma(z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} - {\frac{\omega - \omega_{0}}{v_{g}}z}} \right\rbrack} \right\}\ {\mathbb{d}\omega}}}} & (B) \\ {D = {{- \frac{2\;\pi\; c}{\lambda^{2}}}\beta^{''}}} & (C) \end{matrix}$

-   S(ω): spectral density function of known wide-band spectral light -   d: relative delay difference between first and second interference     signal detecting means 19 and 23 -   ω₀: central frequency of wide-band spectral light -   ν_(g): group velocity of the measured optical fiber 11 (can be     calculated from refractive index or the like of the measured optical     fiber 11) -   z₀: length of measuring section extending from the incident end 12     of the measured optical fiber 11 -   λ: central wavelength of wide-band spectral light -   K: constant -   |I₁|², |I₂|²: squares of absolute values of I₁ and I₂     |I²|:average of |I₁|² and |I₂|²  [Numerical Expression 113]

The wide-band spectral light generating means 14 is realized by, for example, a super-luminescent diode and generates a wide-band spectral continuous light having a spectral width of about several ten nanometers. The variable bandpass filter means 15 cuts out spectral light having a central frequency ω₀ and a bandwidth Δ_(ω). The central frequency ω₀ of the spectral light can be varied by the data acquisition control means 26. It is assumed that the spectral density function S(ω) of the continuous light cut out by the variable bandpass filter means 15 is known in advance. In general, an output from the variable bandpass filter means 15 can be approximated by a Gaussian profile and represented as the following Numerical Expression (114).

$\begin{matrix} {{S(\omega)} = \frac{\exp\left\lbrack \frac{- \left( {\omega - \omega_{0}} \right)^{2}}{2\left( {\Delta\;\omega} \right)^{2}} \right\rbrack}{\sqrt{2\;\pi}\Delta\;\omega}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 114} \right\rbrack \end{matrix}$ Here, in Numerical Expression (114), reference symbol ω₀ is a central frequency which can be controlled by the data acquisition control means 26. Reference symbol Δ_(ω) is a bandwidth.

A signal being proportional to a Rayleigh-scattered light amplitude generated at a first halfway point in the measuring section of the measured optical fiber 11 by the wide-band spectral light (reference light 90 b) being incident on the measured optical fiber 11 and a signal being proportional to a Rayleigh-scattered light amplitude generated at a second halfway point included in the measuring section of the measured optical fiber 11 and being different from the first halfway point are observed. In order to realize this, in Embodiment 2-1, the Rayleigh-scattered light 91 a generated by the wide-band spectral light (reference light 90 b) being incident on the measured optical fiber 11 and the wide-band spectral light (reference light 90 c) transmitted through the delay means 20 having delay time equal to propagation time corresponding to twice the measuring section and the variable delay means 27 are guided to the first interference signal detecting means 19 and the second interference signal detecting means 23 which have a relative optical path length being equal to a distance between the first halfway point and the second halfway point to detect signals as the interference signal I₁ being proportional to the Rayleigh-scattered light amplitude from the first halfway point and the interference signal I₂ being proportional to the Rayleigh-scattered light amplitude from the second halfway point.

More specifically, after the wide-band spectral light (reference light 90 a) passes through the variable bandpass filter means 15, the wide-band spectral light (reference light 90 a) is split by two by means of the optical splitting means 16. One of the split lights is guided to the measured optical fiber 11. At this time, the Rayleigh-scattered light 91 a caused by a fluctuation in refractive index of glass at an arbitrary point of the measured optical fiber 11 is generated, and the generated Rayleigh-scattered light 91 a is guided to the first and second interference signal detecting means 19 and 23 through the optical circulate means 17. The other wide-band spectral light (reference light 90 c) of the lights split by the optical splitting means 16 is guided to the first and second interference signal detecting means 19 and 23 through the delay means 20 having delay time equal to propagation time corresponding to twice the measuring section length and the variable delay means 27. The delay means 20 having delay time equal to the propagation time corresponding to twice the measuring section length is set so as to make the optical path lengths of the following optical paths 1 and 2 equal to each other.

Optical path 1: optical splitting means 16→optical circulate means 17→arbitrary reflection point (of Rayleigh scattering) of measured optical fiber 11→optical circulate means 17→first and second interference signal detecting means 19 and 23. Optical path 2: optical splitting means 16→first and second interference signal detecting means 19 and 23.

More specifically, this can be realized by changing connections of a fixed delay line of an optical fiber. The length of the optical fiber is changed to make it possible to arbitrarily change a measuring section.

The reference light 90 c and the Rayleigh-scattered light 91 a are split by two by the optical splitting means 21 and 18, respectively. The first interference signal detecting means 19 and the second interference signal detecting means 23 which have a relative optical path length difference equal to a distance between the first halfway point and the second halfway point of the measured optical fiber 11 detect interference signals generated by the reference lights 90 c and 90 d and the Rayleigh-scattered lights 91 a and 91 b. When the delay means 20 having delay time equal to propagation time corresponding to twice the measuring section length and the variable delay means 27 are inserted into the path of the reference light 90 c, the interference signal I₁ being proportional to the Rayleigh-scattered light amplitude from the first halfway point of the measured optical fiber 11 is output from the first interference signal detecting means 19. On the other hand, from the second interference signal detecting means 23, the interference signal I₂ being proportional to the Rayleigh-scattered light amplitude is output from the second halfway point of the measured optical fiber 11.

In this embodiment, the interference signals detected by the first interference signal detecting means 19 and the second interference signal detecting means 23 are finally guided to the absolute-value-square calculating means 51, the data acquiring/storing means 24, and the data calculating means to make it possible to derive a wavelength dispersion of a predetermined section by a theoretical expression given by a knowledge (will be described later).

The first and second interference signal detecting means 19 and 23 can be constituted by a heterodyne detecting technique or a homodyne detecting technique which is a known technique. As constituent elements thereof, several variations can be conceived.

The reason why the above relational expression (Numerical Expression (112)) is established will be described below. The waveform of the wide-band spectral light can be represented as the following Numerical Expression (115). Ψ(t)=∫_(−∞) ^(∞) A(ω)exp{i[k(ω)z−ωt]}dω  [Numerical Expression 115]

In this expression, reference symbol A(ω) denotes an electric field amplitude component which is Fourier-dissolved. The following Numerical Expression (116) is established.

$\begin{matrix} {{S(\omega)} = \frac{{{A(\omega)}}^{2}}{2\;\pi}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 116} \right\rbrack \end{matrix}$ In Numerical Expression (116), reference symbol k is a propagation constant, and it is assumed that the path of the reference light 90 c is nondispersive. At this time, the following Numerical Expression (117) is established.

$\begin{matrix} {{k(\omega)} = {\frac{\omega_{0}}{v_{p}} + \frac{\omega - \omega_{0}}{v_{g}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 117} \right\rbrack \end{matrix}$

The Rayleigh-scattered lights 91 a from respective points z_(i) of the optical fiber is expressed by the following Numerical Expression (118) in the inputs to the first and second interference signal detecting means 19 and 23: Ψ_(r) (t)=α(z _(i))·∫_(−∞) ^(∞) A(ω)exp{i[2β(ω)z _(i) −ωt]}dω  [Numerical Expression 118] In this equation, a(z_(i)) denotes a scattering coefficient (local refractive index) in z_(i). In Numerical Expression (118), β(ω) denotes a propagation constant of the optical fiber. According to the Taylor expansion, β(ω) can be represented as the following Numerical Expression (119):

$\begin{matrix} {{\beta(\omega)} = {\frac{\omega_{0}}{v_{p}} + \frac{\omega - \omega_{0}}{v_{g}} + {\frac{\beta^{''}}{2}\left( {\omega - \omega_{0}} \right)^{2}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 119} \right\rbrack \end{matrix}$ In Numerical Expression (119), reference symbol v_(p) denotes a phase velocity, reference numeral ν_(g) denotes a group velocity of light, and a value β″ has a relationship given by the following Numerical Expression (120) to the D. When the value β″ is calculated, the D can be calculated:

$\begin{matrix} {D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 120} \right\rbrack \end{matrix}$

When a total sum of Numerical Expression (118) with respect to z_(i) is calculated, an amplitude of the Rayleigh scattering is obtained as Numerical Expression (121):

$\begin{matrix} {{\Psi_{total}(t)} = {\sum\limits_{i}{{a\left( z_{i} \right)} \cdot {\int_{- \infty}^{\infty}{{A(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{2{\beta(\omega)}z_{i}} -} \\ {\omega\; t} \end{bmatrix}} \right\}\ {\mathbb{d}\omega}}}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 121} \right\rbrack \end{matrix}$

Outputs from the first interference signal detecting means 19 and the second interference signal detecting means 23 are given by the following Numerical Expression (122): I ₁=

Ψ_(total)(t)+Ψ(t)|²

I ₂=

Ψ_(total) (t)+Ψ(t+d/c)|²

  [Numerical Expression 122] (d expresses length of reference optical path, c expresses velocity of light, and < > expresses time averaging operation.). For this reason, when Numerical Expression (115) and Numerical Expression (121) are assigned to Numerical Expression (122) to neglect a DC component, the following Numerical Expression 123 is given:

$\begin{matrix} \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 123} \right\rbrack & \; \\ \begin{matrix} {I_{1} = {{\begin{matrix} {{\sum\limits_{i}{a{\left( z_{i} \right) \cdot {\int_{- \infty}^{\infty}{{A(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack {{2{\beta(\omega)}z_{i}} - {\omega\; t}} \right\rbrack} \right\}{\mathbb{d}\omega}}}}}} +} \\ {\int_{- \infty}^{\infty}{{A(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack {{{k\left( \overset{.}{\omega} \right)}z_{0}} - {\omega\; t}} \right\rbrack} \right\}{\mathbb{d}\omega}}} \end{matrix}}^{2}}} \\ {= {\sum\limits_{i}{{a\left( z_{i} \right)}{\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}}{2}\left( {\omega - \omega_{0}} \right)^{2}z_{i}} +} \\ {{\frac{\omega - \omega_{0}}{v_{g}}z_{i}} + {\frac{\omega_{0}}{v_{p}}z_{i}}} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}}} \\ {= {\sum\limits_{i}{{a\left( z_{i} \right)}{\Gamma\left( z_{i} \right)}}}} \end{matrix} & (A) \\ {\mspace{79mu}\begin{matrix} {I_{2} = {{\begin{matrix} {{\sum\limits_{i}{a{\left( z_{i} \right) \cdot {\int_{- \infty}^{\infty}{{A(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack {{2{\beta(\omega)}z_{i}} - {\omega\; t}} \right\rbrack} \right\}{\mathbb{d}\omega}}}}}} +} \\ {\int_{- \infty}^{\infty}{{A(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack {{{k(\omega)}\left( {z + d} \right)} - {\omega\; t}} \right\rbrack} \right\}{\mathbb{d}\omega}}} \end{matrix}}^{2}}} \\ {= {\sum\limits_{i}{{a\left( z_{i} \right)}{\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}}{2}\left( {\omega - \omega_{0}} \right)^{2}\left( {z_{i} - d} \right)} +} \\ {{\frac{\omega - \omega_{0}}{v_{g}}\left( {z_{i} - d} \right)} +} \\ {\frac{\omega_{0}}{v_{p}}\left( {z_{i} - d} \right)} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}}} \\ {= {\sum\limits_{i}{{a\left( z_{i} \right)}{\Gamma\left( {z_{i} - d} \right)}}}} \end{matrix}} & (B) \\ {\mspace{79mu}{where}} & \; \\ {\mspace{79mu}{{\Gamma(z)} = {{\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}}z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{S(\omega)\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}}{2}\left( {\omega - \omega_{0}} \right)^{2}z} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}}} & \begin{bmatrix} {Numerical} \\ {Expression} \\ 124 \end{bmatrix} \end{matrix}$ In general, Γ(z) is called a coherence function of wide-band spectral light. In this case, when it is assumed that β″z/2 changes by only quantity which is considerably smaller than 2π with reference to a change of z which is almost equal to d, Γ(z) may be expressed by the following numerical expression:

$\begin{matrix} {{\Gamma(z)} \cong {{\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}}z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{S(\omega)\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}} & \begin{bmatrix} {Numerical} \\ {Expression} \\ 125 \end{bmatrix} \end{matrix}$

A real part and an imaginary part of (A) and (B) of Numerical Expression (123) are expressed as |Γ(z)|α_(i) cos θ_(i) and |Γ(z)|α_(i) sin θ_(i). In this case, a vector (α_(i) cos θ_(i), α_(i) sin θ_(i)) represents a random phaser, and α_(i) represents random variables each expressing a length of the phaser. It is assumed that an average of random variables is given by: α  [Numerical Expression 126] and that a second moment is given by: ( α)²  [Numerical Expression 127] Furthermore, θ_(i) represents random variables each expressing a phase of the phaser. The random variables are uniformly distributed within the range of −π to π. The variables α_(i) and θ_(i) are naturally independent of each other. When the expressions of the real part and the imaginary part are used, the following expressions are obtained:

$\begin{matrix} \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 128} \right\rbrack & \; \\ \begin{matrix} {{I_{1}}^{2} = {\left\{ {\sum\limits_{i}{{{\Gamma\left( z_{i} \right)}}\alpha_{i}\cos\;\theta_{i}}} \right\}^{2} + \left\{ {\sum\limits_{i}{{{\Gamma\left( z_{i} \right)}}\alpha_{i}\sin\;\theta_{i}}} \right\}^{2}}} \\ {= {{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}\alpha_{i}\alpha_{j}\cos\;\theta_{i}\cos\;\theta_{j}}}} +}} \\ {\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}\alpha_{i}\alpha_{j}\sin\;\theta_{i}\sin\;\theta_{j}}}} \end{matrix} & (A) \\ \begin{matrix} {{I_{2}}^{2} = {\left\{ {\sum\limits_{i}{{{\Gamma\left( {z_{i} - d} \right)}}\alpha_{i}\cos\;\theta_{i}}} \right\}^{2} + \left\{ {\sum\limits_{i}{{{\Gamma\left( {z_{i} - d} \right)}}\alpha_{i}\sin\;\theta_{i}}} \right\}^{2}}} \\ {= {{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}\alpha_{i}\alpha_{j}\cos\;\theta_{i}\cos\;\theta_{j}}}} +}} \\ {\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}\alpha_{i}\alpha_{j}\sin\;\theta_{i}\sin\;\theta_{j}}}} \end{matrix} & (B) \end{matrix}$ Therefore, the following expression is obtained:

$\begin{matrix} {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} = {{\sum\limits_{i}{\sum\limits_{j}{\sum\limits_{k}{\sum\limits_{n}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{k} - d} \right)}}{{\Gamma\left( {z_{n} - d} \right)}}{\overset{\_}{\alpha_{i}\alpha_{j}\alpha_{k}\alpha_{n}} \cdot \overset{\_}{\cos\;\theta_{i}\cos\;\theta_{j}\cos\;\theta_{k}\cos\;\theta_{n}}}}}}}} + {\sum\limits_{i}{\sum\limits_{j}{\sum\limits_{k}{\sum\limits_{n}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{k} - d} \right)}}{{\Gamma\left( {z_{n} - d} \right)}}{\overset{\_}{\alpha_{i}\alpha_{j}\alpha_{k}\alpha_{n}} \cdot \overset{\_}{\cos\;\theta_{i}\cos\;\theta_{j}\sin\;\theta_{k}\sin\;\theta_{n}}}}}}}} + {\sum\limits_{i}{\sum\limits_{j}{\sum\limits_{k}{\sum\limits_{n}{{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}{{\Gamma\left( z_{k} \right)}}{{\Gamma\left( z_{n} \right)}}{\overset{\_}{\alpha_{i}\alpha_{j}\alpha_{k}\alpha_{n}} \cdot \overset{\_}{\cos\;\theta_{i}\cos\;\theta_{j}\sin\;\theta_{k}\sin\;\theta_{n}}}}}}}} + {\sum\limits_{i}{\sum\limits_{j}{\sum\limits_{k}{\sum\limits_{n}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{k} - d} \right)}}{{\Gamma\left( {z_{n} - d} \right)}}{\overset{\_}{\alpha_{i}\alpha_{j}\alpha_{k}\alpha_{n}} \cdot \overset{\_}{\sin\;\theta_{i}\sin\;\theta_{j}\sin\;\theta_{k}\sin\;\theta_{n}}}}}}}}}} & \begin{bmatrix} {Numerical} \\ {Expression} \\ 129 \end{bmatrix} \end{matrix}$ Since the variables θ_(i) are uniformly distributed within the range of −π to π, a term which is not zero in the right-hand side of Numerical Expression (129) appears in only the following case: 1. In case of i=j=k=n, cos θ_(i) cos θ_(j) cos θ_(k) cos θ_(n) = cos⁴θ_(i) =⅜ cos θ_(i) cos θ_(j) cos θ_(k) cos θ_(n) = cos²θ_(i) sin²θ_(i) =⅛ sin θ_(i) sin θ_(j) sin θ_(k) sin θ_(n) = sin⁴θ_(i) =⅜  [Numerical Expression 130] 2. In case of i=j and k=n, cos θ_(i) cos θ_(j) cos θ_(k) cos θ_(n) = cos²θ_(i) · cos²θ_(k) =¼ cos θ_(i) cos θ_(j) sin θ_(k) sin θ_(n) = cos²θ_(i) · sin²θ_(k) =¼ sin θ_(i) sin θ_(j) sin θ_(k) sin θ_(n) = sin²θ_(i) ·sin²θ_(k)¼ 3. In case of i=k and j=n, cos θ_(i) cos θ_(j) cos θ_(k) cos θ_(n) = cos²θ_(i) · cos²θ_(j) =¼ cos θ_(i) cos θ_(j) cos θ_(k) cos θ_(n) = cos θ_(i) sin θ_(i) · cos θ_(j) sin θ_(j) =¼ sin θ_(i) sin θ_(j) sin θ_(k) sin θ_(n) = sin²θ_(i) · sin²θ_(j) =¼ 4. In case of i=n, j=k, cos θ_(i) cos θ_(j) cos θ_(k) cos θ_(n) = cos²θ_(i) · cos²θ_(j) =¼ cos θ_(i) cos θ_(j) sin θ_(k) sin θ_(n) = cos θ_(i) sin θ_(i) · cos θ_(j) sin θ_(j) =¼ sin θ_(i) sin θ_(j) sin θ_(k) sin θ_(n) = sin²θ_(i) · sin²θ_(j) =¼ In all other cases except for the above cases, the following expression is obtained: cos θ_(i) cos θ_(j) cos θ_(k) cos θ_(n) =0 cos θ_(i) cos θ_(j) sin θ_(k) sin θ_(n) =0 sin θ_(i) sin θ_(j) sin θ_(k) sin θ_(n) =0  [Numerical Expression 131]

Although the case of i=j and k=n, the case of i=k and j=n, and the case of i=n and j=k appear at a equally probability (N(N−1) times), the case of i=j=k=n appears at a low probability (N times). For this reason, the case can be neglected (becomes zero when Σ is replaced with integration). For this reason, Numerical Expression (129) can be expressed as:

$\begin{matrix} \begin{matrix} {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} = {\frac{1}{4}{\sum\limits_{i}{\sum\limits_{k}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{k} - d} \right)}}}}}}} \\ {{{{\Gamma\left( {z_{k} - d} \right)}}\overset{\_}{\alpha_{i}\alpha_{i}\alpha_{k}\alpha_{k}}} +} \\ {\frac{1}{4}{\sum\limits_{i}{\sum\limits_{k}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{k} - d} \right)}}}}}} \\ {{{{\Gamma\left( {z_{k} - d} \right)}}\overset{\_}{\alpha_{i}\alpha_{i}\alpha_{k}\alpha_{k}}} +} \\ {\frac{1}{4}{\sum\limits_{i}{\sum\limits_{k}{{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( z_{k} \right)}}}}}} \\ {{{{\Gamma\left( z_{k} \right)}}\overset{\_}{\alpha_{i}\alpha_{i}\alpha_{k}\alpha_{k}}} +} \\ {\frac{1}{4}{\sum\limits_{i}{\sum\limits_{k}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{k} - d} \right)}}}}}} \\ {{{{\Gamma\left( {z_{k} - d} \right)}}\overset{\_}{\alpha_{i}\alpha_{i}\alpha_{k}\alpha_{k}}} +} \\ {\frac{1}{4}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{i} - d} \right)}}}}}} \\ {{{{\Gamma\left( {z_{j} - d} \right)}}\overset{\_}{\alpha_{i}\alpha_{j}\alpha_{i}\alpha_{j}}} +} \\ {\frac{1}{4}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}} \\ {{{{\Gamma\left( {z_{j} - d} \right)}}\overset{\_}{\alpha_{i}\alpha_{j}\alpha_{j}\alpha_{j}}} +} \\ {\frac{1}{4}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}} \\ {{{{\Gamma\left( {z_{i} - d} \right)}}\overset{\_}{\alpha_{i}\alpha_{j}\alpha_{j}\alpha_{i}}} +} \\ {\frac{1}{4}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}} \\ {{{\Gamma\left( {z_{i} - d} \right)}}\overset{\_}{\alpha_{i}\alpha_{j}\alpha_{j}\alpha_{i}}} \\ {= {{\frac{\left( \overset{\_}{\alpha^{2}} \right)^{2}}{4}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}} +}} \\ {{\frac{\left( \overset{\_}{\alpha^{2}} \right)^{2}}{4}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}} +} \\ {{\frac{\left( \overset{\_}{\alpha^{2}} \right)^{2}}{4}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}} +} \\ {{\frac{\left( \overset{\_}{\alpha^{2}} \right)^{2}}{4}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}} +} \\ {{\frac{\left( \overset{\_}{\alpha^{2}} \right)^{2}}{4}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}} +} \\ {{\frac{\left( \overset{\_}{\alpha^{2}} \right)^{2}}{4}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}} +} \\ {{\frac{\left( \overset{\_}{\alpha^{2}} \right)^{2}}{4}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}} +} \\ {\frac{\left( \overset{\_}{\alpha^{2}} \right)^{2}}{4}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}} \end{matrix} & \begin{bmatrix} {Numerical} \\ {Expression} \\ 132 \end{bmatrix} \end{matrix}$ Four terms of the first half of Numerical Expression (132) are equal to each other, and four terms of the second half are equal to each other. Since a delay d in the four terms of the first half acts on only z_(j), the delay d can be removed without influence. Therefore, the following expression is established:

$\begin{matrix} {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} = {{\left( \overset{\_}{\alpha^{2}} \right)^{2}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( z_{j} \right)}}}}}} + {\left( \overset{\_}{\alpha^{2}} \right)^{2}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}}}} & \begin{bmatrix} {Numerical} \\ {Expression} \\ 133 \end{bmatrix} \end{matrix}$

From (A) and (B) of Numerical Expression (128), the following Numerical Expression (134) is established:

$\begin{matrix} {\begin{matrix} {\overset{\_}{{I_{1}}^{2}} = {{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}\overset{\_}{\alpha_{i}\alpha_{j}\cos\;\theta_{i}\cos\;\theta_{j}}}}} +}} \\ {\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}\overset{\_}{\alpha_{i}\alpha_{j}\sin\;\theta_{i}\sin\;\theta_{j}}}}} \\ {= {\overset{\_}{\alpha^{2}}{\sum\limits_{i}{{\Gamma\left( z_{i} \right)}}^{2}}}} \\ {\equiv \overset{\_}{{I}^{2}}} \end{matrix}\begin{matrix} {\overset{\_}{{I_{2}}^{2}} = {{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}\overset{\_}{\alpha_{i}\alpha_{j}\cos\;\theta_{i}\cos\;\theta_{j}}}}} +}} \\ {\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}\overset{\_}{\alpha_{i}\alpha_{j}\sin\;\theta_{i}\sin\;\theta_{j}}}}} \\ {= {\overset{\_}{\alpha^{2}}{\sum\limits_{i}{{\Gamma\left( z_{i} \right)}}^{2}}}} \\ {\equiv \overset{\_}{{I}^{2}}} \end{matrix}} & \begin{bmatrix} {Numerical} \\ {Expression} \\ 134 \end{bmatrix} \end{matrix}$ Therefore, the following numerical expression is established:

$\begin{matrix} {\overset{\_}{{I_{1}}^{2}} = {{+ \left( \overset{\_}{\alpha^{2}} \right)^{2}}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( z_{j} \right)}}}}}}} & \begin{bmatrix} {Numerical} \\ {Expression} \\ 135 \end{bmatrix} \end{matrix}$ Since this is equal to the first term of the right-hand side of Numerical Expression (133), according to Numerical Expression (133), the following numerical expression is obtained:

$\begin{matrix} {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} = {\left( \overset{\_}{{I}^{2}} \right)^{2} + {\left( \overset{\_}{\alpha^{2}} \right)^{2}{\sum\limits_{i}{\sum\limits_{j}{{{\Gamma\left( z_{i} \right)}}{{\Gamma\left( {z_{i} - d} \right)}}{{\Gamma\left( z_{j} \right)}}{{\Gamma\left( {z_{j} - d} \right)}}}}}}}} & \begin{bmatrix} {Numerical} \\ {Expression} \\ 136 \end{bmatrix} \end{matrix}$ When Σ is replaced with integration, |I ₁ | ² |I ₂|² =( |I| ² )²+( α² )²∫_(−∞) ^(∞)|Γ(z _(i))∥Γ(z _(i) −d)∥Γ(z _(j))∥Γ(z _(j) −d)|dz _(i) dz=( |I| ² )²+( α² )²{∫_(−∞) ^(∞)|Γ(z)|·|Γ(z−d)|dz} ²  [Numerical Expression 137] is established. In this case, when (α²) ²  [Numerical Expression 138] is replaced with K, Numerical Expression 137 is equal to (A) of Numerical Expression (112). In this case, Γ is given by:

$\begin{matrix} {{\Gamma(z)} \cong {{\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}}z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}} & \begin{bmatrix} {Numerical} \\ {Expression} \\ 139 \end{bmatrix} \end{matrix}$ However, since all Γ of (A) of Numerical Expression (112) calculate absolute values, contribution of the following numerical expression:

$\begin{matrix} {\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}}z} \right\rbrack} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 140} \right\rbrack \end{matrix}$ is eliminated, Numerical Expression 139 may be expressed as Numerical Expression 141:

$\begin{matrix} {{\Gamma(z)} = {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 141} \right\rbrack \end{matrix}$ Therefore, it was proved that the following Numerical Expression (142) is established:

$\begin{matrix} \begin{matrix} {\overset{\_}{\left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right)\left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right)} = {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} -}} \\ {\overset{\_}{{I}^{2}}} \\ {= {K\begin{Bmatrix} {\int_{- \infty}^{\infty}{{{\Gamma(z)}} \cdot}} \\ {{{\Gamma\left( {z - d} \right)}}{\mathbb{d}z}} \end{Bmatrix}^{2}}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 142} \right\rbrack \\ {\;{{where},}} & \; \\ {{\Gamma(z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}} & \; \end{matrix}$

(Embodiment 2-2) FIG. 5 is an explanatory block diagram showing a wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 2-2 of the present invention. This embodiment will describe a case in which an optical waveguide is a measured optical fiber 11. The optical waveguide is not limited to the measured optical fiber 11, and a slab type or channel type optical waveguide which is used in a planar lightwave circuit and which confines light in a predetermined medium to propagate the light can also be applied. The same reference numerals as in Embodiment 2-1 denote the same parts in Embodiment 2-2.

As shown in FIG. 5, a wavelength dispersion of a section (to be referred as a measuring section hereinafter in the embodiment) extending from an incident end 12 of a measured optical fiber 11 to an arbitrary halfway point 13 is measured. Wide-band spectral light (reference light 90 a) having a known spectral density function S(ω) from a wide-band spectral light generating means 14 is incident on the measured optical fiber 11 through a variable bandpass filter means 15, an optical splitting means 16, and an optical circulate means 17. Rayleigh-scattered light 91 a generated at a halfway point of the measured optical fiber 11 by the wide-band spectral light (reference light 90 b) being incident on the measured optical fiber 11 is input to an interference signal detecting means 30 serving as a first interference signal detecting means and a second interference signal detecting means through the optical circulate means 17. Wide-band spectral light (reference light 90 c) split by the optical splitting means 16 is input to the interference signal detecting means 30 through a delay means 20 having delay time equal to propagation time corresponding to twice a length of a measuring section from the incident end 12 of the measured optical fiber 11 and a variable delay means 27 interlocked with a data acquisition control means 26.

In the interference signal detecting means 30, signals being proportional to Rayleigh-scattered light amplitudes are detected to be input to an absolute-value-square calculating means 51, a data acquiring/storing means 24, and a data calculating means 25. In the absolute-value-square calculating means 51, the data acquiring/storing means 24, and the data calculating means 25, a wavelength dispersion of the section extending from the incident end 12 of the measured optical fiber 11 to the arbitrary halfway point 13 is calculated by a correlational function of signals being proportional to the Rayleigh-scattered light intensities generated at a plurality of halfway points. The variable delay means 27 and the data acquiring/storing means 24 are controlled by the data acquisition control means 26. The data acquisition control means 26 records outputs from the interference signal detecting means 30 as a function of delay time τ(=d/c) on the data acquiring/storing means 24 while changing quantity of delay of the variable delay means 27. The interference signals detected by the interference signal detecting means 30 are stored by the data acquiring/storing means 24. In the absolute-value-square calculating means 51, a square of an absolute value is calculated from I(τ) recorded as a function of the delay time τ. In the data calculating means 25, a plurality of numerical values at two arbitrary points time d/c distant from each other are extracted as |I₁|² and |I₂|² from |I(τ)|² which is a square of the obtained absolute value according to a condition given by the following Numerical Expression (143), and, by using relational expressions ((A), (B), and (C) of Numerical Expression 112) in the Embodiment 2-1, a dispersion of the measured optical fiber 11 in the predetermined section is calculated:

$\begin{matrix} \begin{matrix} {{{\overset{\_}{I}}_{1}}^{2} = {{I\left( \tau_{i} \right)}}^{2}} \\ {{{\overset{\_}{I}}_{2}}^{2} = {{I\left( {\tau_{i} + \frac{d}{c}} \right)}}^{2}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 143} \right\rbrack \end{matrix}$

The wavelength dispersion measuring method for an optical fiber according to the embodiments is executed by a wavelength dispersion measuring apparatus for an optical fiber including the same functions as those of, more specifically, a computer such as a personal computer on the basis of a predetermined measuring program in advance. The measuring program can be recorded on a recording medium such as a CD which can be read by the wavelength dispersion measuring apparatus for an optical fiber.

More specifically, the measuring program is a measuring program used in a wavelength dispersion measuring apparatus for an optical fiber which can measure a wavelength dispersion of a measuring section of an optical fiber by using a scattering phenomenon in which a scattering coefficient at a local point on the optical fiber is independent of time, and is to cause the wavelength dispersion measuring apparatus for an optical fiber according to the embodiments of the present invention to receive an interference signal I₁ being proportional to a scattered light amplitude generated by wide-band spectral light being incident on the optical fiber at a first halfway point and an interference signal I₂ being proportional to a scattered light amplitude generated at a second halfway point different from the first halfway point as numerical values and to calculate a correlational function of squares of absolute values of the interference signals I₁ and I₂, thereby to calculate a wavelength dispersion D in the measuring section by theoretical expressions ((A), (B), and (C) of Numerical Expression 112) in Embodiment 2-1.

In the wavelength dispersion measuring apparatus for an optical fiber according to Embodiment 2-1 of the present invention, when the first interference signal and the second interference signal are represented by I₁ and I₂, respectively, a correlational function between the first interference signal and the second interference signal detected by the first interference signal detecting means 19 and the second interference signal detecting means 23 is given by the following Numerical Expression (144) expressing an ensemble average of products thereof: I₁I₂*  [Numerical Expression 144] where ⋆ is a complex conjugate sign.

According to the theory, a correlational function between the first interference signal I₁ and the second interference signal I₂ can be associated with each other by the dispersion D in the measuring section and the following relational expression (Numerical Expression (145)):

$\begin{matrix} \; & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 145} \right\rbrack \\ {\begin{matrix} {\frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{1}}^{2}}} = \frac{\overset{\_}{I_{1}I_{2}^{*}}}{\overset{\_}{{I_{2}}^{2}}}} \\ {= \frac{\int_{- \infty}^{\infty}{{{\Gamma(Z)}}{\Gamma^{*}\left( {Z - d} \right)}\ {\mathbb{d}Z}}}{\int_{- \infty}^{\infty}{{{\Gamma(Z)}}^{2}\ {\mathbb{d}Z}}}} \end{matrix}{{where},}} & (A) \\ {{\Gamma(z)} \cong {{\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}}z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}} & (B) \\ {D = {{- \frac{2\;\pi\; c}{\lambda^{2}}}\beta^{''}}} & (C) \end{matrix}$

The above-mentioned wavelength dispersion measuring method for an optical fiber is also executed by a wavelength dispersion measuring apparatus for an optical fiber including the same functions as those of, more specifically, a computer such as a personal computer on the basis of a predetermined measuring program in advance. The measuring program can be recorded on a recording medium such as a CD which can be read by the wavelength dispersion measuring apparatus for an optical fiber.

More specifically, the measuring program is a measuring program used in a wavelength dispersion measuring apparatus for an optical fiber which can measure a wavelength dispersion of a measuring section of an optical fiber by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, and is to cause the wavelength dispersion measuring apparatus for an optical fiber to receive an interference signal I₁ being proportional to a scattered light amplitude generated by wide-band spectral light being incident on the optical fiber at a first halfway point and an interference signal I₂ being proportional to a scattered light amplitude generated at a second halfway point different from the first halfway point as numerical values and to calculate a correlational function of the interference signals I₁ and I₂, thereby to calculate a wavelength dispersion D in the measuring section by (A), (B), and (C) of the above Numerical Expression 145.

The present invention is not directly limited to the above embodiment. In an execution phase, the invention can be embedded by modifying constituent elements within departing from the spirit and scope of the invention. Appropriate combinations of a plurality of constituent elements disclosed in the above embodiment can constitute various inventions. For example, several constituent elements may be omitted from all the constituent elements disclosed in the embodiments. Furthermore, the constituent elements across different embodiments may be appropriately combined to each other.

Embodiment 3

FIG. 6 shows a configuration of an embodiment of the present invention. This embodiment will describe a case in which an optical waveguide is a measured optical fiber 11. The optical waveguide is not limited to the measured optical fiber 11, and a slab type or channel type optical waveguide which is used in a planar lightwave circuit and which confines light in a predetermined medium to propagate the light can also be applied. The same reference numerals as in the other embodiments denote the same constituent elements in this embodiment.

In the embodiment, in order to measure Rayleigh-scattered light 91 a from an optical fiber, an optical frequency domain reflectometry method (to be referred to as an OFDR method hereinafter) is used. It is a known technique (for example, see Patent Document 1) that an amplitude or an intensity of the Rayleigh-scattered light 91 a from the measured optical fiber 11 can be measured as a function of location by the OFDR method.

Patent Document 1: Japanese Patent No. 2907350

In the embodiment, as a means for observing the Rayleigh-scattered light 91 a, a low-coherence reflecting method explained in another embodiment can also be used. The measuring means is not limited to a specific means. In the OFDR method, a narrow-line-width laser is used as a light source serving as a narrow-line-width frequency modulated light generating means 61.

In this case, FIG. 7 shows a schematic diagram showing a frequency modulation waveform in the OFDR method, and FIG. 8 shows a schematic diagram showing a spectral density function in the OFDR method.

In the narrow-line-width frequency modulated light generating means 61, the oscillation frequency thereof is linearly changed with respect to time as shown in FIG. 7. A spectral density function S(ω) at this time has a rectangular shape as shown in FIG. 8.

As shown in FIG. 6, a wavelength dispersion of a section (to be referred to as a measuring section hereinafter in the embodiment) extending from an incident end 12 of the measured optical fiber 11 to a halfway point 13 is measured. Wide-band spectral light (reference light 90 a) having a known spectral density function S(ω) from the narrow-line-width frequency modulated light generating means 61 is incident on the measured optical fiber 11 through an optical splitting means 16, and an optical circulate means 17. Rayleigh-scattered light 91 a generated at a halfway point of the measured optical fiber 11 by the wide-band spectral light (reference light 90 b) being incident on the measured optical fiber 11 is input to an interference signal detecting means 30 serving as a first interference signal detecting means and a second interference signal detecting means through the optical circulate means 17. Wide-band spectral light (reference light 90 c) split by the optical splitting means 16 is input to the interference signal detecting means 30.

In the interference signal detecting means 30, signals being proportional to Rayleigh-scattered light amplitudes are detected to be input to a data acquiring/storing means 24, an every-distance Rayleigh-scattered light amplitude analyzing means 101, an every-distance Rayleigh-scattered light amplitude data storing means 102, a absolute-value-square calculating means 51, and a data calculating means 25. The every-distance Rayleigh-scattered light amplitude analyzing means 101 analyzes the interference signals detected by the interference signal detecting means 30 on the basis of the principle of OFDR to calculate Rayleigh-scattered light amplitudes at the respective points of the measured optical fiber 11 and record the Rayleigh-scattered light amplitudes as functions of delay time from the incident end 12 of the measured optical fiber 11 to the respective points on the every-distance Rayleigh-scattered light amplitude data storing means 102. In the absolute-value-square calculating means 51 and the data calculating means 25, a wavelength dispersion of the section extending from the incident end 12 of the measured optical fiber 11 to the arbitrary halfway point 13 is calculated by a correlational function of signals being proportional to Rayleigh-scattered light intensities generated at a plurality of halfway points. The narrow-line-width frequency modulated light generating means 61 and the data acquiring/storing means 24 are controlled by the data acquisition control means 26. The data acquisition control means 26 records outputs from the interference signal detecting means 30 on the data acquiring/storing means 24 at constant intervals of time while controlling a modulation waveform of the narrow-line-width frequency modulated light generating means 61. The interference signals I detected by the interference signal detecting means 30 are stored by the data acquiring/storing means 24. In the absolute-value-square calculating means 51, a square of an absolute value is calculated from I(τ) recorded as a function of the delay time τ from the incident end 12 of the measured optical fiber 11 to respective points. In the data calculating means 25, a plurality of numerical values at two arbitrary points time d/c distant from each other are extracted from |I(τ)|² which is a square of the obtained absolute value according to a condition given by the following Numerical Expression (146), and, by using relational expressions ((A), (B), and (C) of Numerical Expression 148), a dispersion of the measured optical fiber 11 in the predetermined section is calculated:

$\begin{matrix} \begin{matrix} {{{\overset{\_}{I}}_{1}}^{2} = {{I\left( \tau_{i} \right)}}^{2}} \\ {{{\overset{\_}{I}}_{2}}^{2} = {{I\left( {\tau_{i} + \frac{d}{c}} \right)}}^{2}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 146} \right\rbrack \end{matrix}$

The embodiment is characterized in that a reference dispersion giving means 71 in which the following Numerical Expression (147) expressing an accumulated wavelength dispersion is known is inserted on a path of the reference light 90 c from the narrow-line-width frequency modulated light generating means 61 to the interference signal detecting means 30 or a path of the Rayleigh-scattered light 91 a from the measured optical fiber 11: β″_(ref)L_(ref)  [Numerical Expression 147] In the other embodiments described above, an equation which associates a correlational function of scattered light with a wavelength dispersion is an even function to the value of the wavelength dispersion. Therefore, although the absolute value of the wavelength dispersion can be measured, a sign of the wavelength dispersion cannot be known. On the other hand, in the embodiment, it can be shown by the same calculations (Numerical Expressions (115) to (142) in Embodiment 2-1) of another embodiment that a relational expression obtained by inserting the reference dispersion giving means 71 is expressed as the following Numerical Expression (148). The embodiment is different from the other embodiment in only that an accumulated wavelength dispersion is expressed by the above Numerical Expression (147) in Γ(z) in Numerical Expressions (115) to (142) in Embodiment 2-1.

$\begin{matrix} \; & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 148} \right\rbrack \\ {\begin{matrix} {\overset{\_}{\left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right)\left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right)} = {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} -}} \\ {\overset{\_}{{I}^{2}}} \\ {= {K\begin{Bmatrix} {\int_{- \infty}^{\infty}{{{\Gamma(z)}} \cdot}} \\ {{{\Gamma\left( {z - d} \right)}}{\mathbb{d}z}} \end{Bmatrix}^{2}}} \end{matrix}} & (A) \\ {{where},} & \; \\ {{\Gamma(z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} \begin{matrix} \frac{{\beta^{''}z_{0}} - {\beta_{ref}^{''}L_{ref}}}{2} \\ {\left( {\omega - \omega_{0}} \right)^{2} -} \end{matrix} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}} & (B) \\ {D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}} & (C) \end{matrix}$

-   S(ω): spectral density function of known wide-band spectral light -   d: relative delay difference -   ω₀: central frequency of wide-band spectral light -   ν_(g): group velocity of the measured optical fiber 11 (can be     calculated from refractive index or the like of the measured optical     fiber 11) -   z₀: length of measuring section extending from the incident end 12     of the measured optical fiber 11 -   λ: central wavelength of wide-band spectral light -   K: constant -   |I₁|², |I₂|²: squares of absolute values of I₁ and I₂     |I|² :average of |I₁|² and |I₂|²[Numerical Expression 149]

Therefore, according to the embodiment, if the wavelength dispersion of the measured optical fiber 11 is zero, a wavelength dispersion corresponding to a value apparently calculated by Numerical Expression (148) is observed. In this manner, according to the embodiment, since a zero point is shifted, a sign of the wavelength dispersion of the measured optical fiber 11 can be detected as an increase or a decrease in absolute value of the apparently measured wavelength dispersion.

The present invention is not directly limited to the example of the above embodiment. In an execution phase, the invention can be embedded by modifying constituent elements within range without departing from the spirit and scope of the invention. Appropriate combinations of a plurality of constituent elements disclosed in the above embodiment can constitute various inventions. For example, several constituent elements may be omitted from all the constituent elements disclosed in the embodiments. Furthermore, the constituent elements across different embodiments may be appropriately combined to each other. In the embodiment, as shown in FIG. 6, a square of an absolute value is calculated from a continuous interference signal I by the interference signal detecting means 30, and a plurality of numerical values at two arbitrary points time d/c distant from each other are extracted as |I₁|² and |I₂|² from the square of the absolute value according to a condition given by the following Numerical Expression (146). However, for example, in the configuration in which the narrow-line-width frequency modulated light generating means 61 is applied in place of the wide-band spectral light generating means 14 in FIG. 1 and the reference dispersion giving means 71 in FIG. 6 is inserted between the delay means 20 and the variable delay means 27, the interference signals I₁ and I₂ may also be independently detected by the first interference signal detecting means 19 and the second interference signal detecting means 23 to give squares of absolute values of the interference signals as |I₁|² and |I₂|², respectively. 

1. A wavelength dispersion measuring method for an optical waveguide capable of measuring a wavelength dispersion in a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, comprising: a step of causing incident light having a known spectral density function S(ω) to be incident on the optical waveguide; a step of observing, by the incident light being incident on the optical waveguide in the step of causing the incident light to be incident on the optical waveguide, a signal being proportional to a scattered light amplitude generated at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and a signal being proportional to a scattered light amplitude generated at a second halfway point being different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way; a step of calculating by a processor a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by a correlational function between the signal being proportional to the scattered light amplitude at the first halfway point and the signal being proportional to the scattered light amplitude at the second halfway point, the signals being observed in the step of observing the signals which are proportional to the scattered light amplitudes; and a step of using the calculated wavelength dispersion to determine the characteristic of the optical waveguide.
 2. The wavelength dispersion measuring method for an optical waveguide according to claim 1, characterized in that the step of observing the signals being proportional to the scattered light amplitudes is the step of observing of Rayleigh-scattered lights as the scattered lights at the first halfway point and the second halfway point.
 3. The wavelength dispersion measuring method for an optical waveguide according to claim 1 or 2, characterized in that the step of observing the signals being proportional to the scattered light amplitudes is the step of guiding scattered light generated by incident light being incident on the optical waveguide and incident light passing through delay means having delay time equal to propagation time corresponding to twice a length between the incident end of the optical waveguide and the first halfway point to first interference signal detecting means and second interference signal detecting means which have a relative delay difference equal to a distance between the first halfway point and the second halfway point to detect an interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and an interference signal I₂ being proportional to a scattered light amplitude from the second halfway point, and the step of calculating the wavelength dispersion is the step of calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by a correlational function between the interference signal I₁ and the interference signal I₂.
 4. The wavelength dispersion measuring method for an optical waveguide according to claim 3, characterized in that variable delay means set on a path of the scattered light or the incident light is arranged, the step of observing the signals being proportional to the scattered light amplitudes is the step of respectively detecting the plurality of interference signals I₁ and I₂ while causing the variable delay means to change quantity of delay, and the step of calculating the wavelength dispersion is the step of calculating a correlational function between the interference signals I₁ and I₂ by calculating the following numerical expression (1) expressing an ensemble average of products of the interference signals I₁and I₂ detected in the step of observing the signals being proportional to the scattered light amplitudes: I₁I₂*  [Numerical Expression 1].
 5. The wavelength dispersion measuring method for an optical waveguide according to claim 4, characterized in that by using the following numerical expression (4) calculated as the correlational function of the interference signals I₁ and I₂: I₁I₂*  [Numerical Expression 4] a dispersion D in the section extending from the incident end of the optical waveguide to the arbitrary point along the way is calculated by the following numerical expression (5): $\begin{matrix} \begin{matrix} {\frac{\overset{\_}{I_{1}I_{2}^{\star}}}{\overset{\_}{{I_{1}}^{2}}} = \frac{\overset{\_}{I_{1}I_{2}^{\star}}}{\overset{\_}{{I_{2}}^{2}}}} \\ {= \frac{\int_{- \infty}^{\infty}{{\Gamma(Z)}{\Gamma^{\star}\left( {Z - d} \right)}{\mathbb{d}Z}}}{\int_{- \infty}^{\infty}{{{\Gamma(Z)}}^{2}{\mathbb{d}Z}}}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 5} \right\rbrack \\ {{where},} & \; \\ {{\Gamma(z)} \cong {{\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}}z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\begin{bmatrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{bmatrix}} \right\}{\mathbb{d}\omega}}}}} & \; \\ {D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}} & \; \end{matrix}$ S(ω): spectral density function of known incident light d: relative delay difference between first and second interference signal detecting means ω₀: central frequency of incident light ν_(g): group velocity of optical waveguide z₀: length of section extending from incident end of the optical waveguide to arbitrary point along the way λ: central wavelength of incident light |I₁|², |I₂|²: squares of absolute values of I₁ and I₂ β: propagation constant of output fiber.
 6. The wavelength dispersion measuring method for an optical waveguide according to claim 3, characterized in that variable delay means set on a path of the scattered light or the incident light is arranged, the step of observing the signals being proportional to the scattered light amplitudes is the step of detecting the interference signals I₁ and I₂ as interference signals I while causing the variable delay means to continuously change delay time τ, and the step of calculating the wavelength dispersion is the step of recording the interference signals I detected in the step of observing the signals being proportional to the scattered light amplitudes as a function I(τ) of τ and calculating a correlational function between the interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and the interference signal I₂ being proportional to a scattered light amplitude from the second halfway point by calculating the following numerical expression (3) expressing an ensemble average of the interference signals I₁ and I₂: I₁I₂*  [Numerical Expression 3] by the following numerical expression (2): I ₁ I ₂*=

I(τ_(i))I(τ_(i) +d/c)

_(i)  [Numerical Expression 2] (where <>_(i) expresses an averaging operation about τi) by using the function I(τ).
 7. A wavelength dispersion measuring apparatus for an optical waveguide capable of measuring a wavelength dispersion in a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, comprising: means for causing incident light having a known spectral density function S(ω) to be incident on the optical waveguide; means for observing, by the incident light from the means for causing the incident light on the optical waveguide, a signal being proportional to a scattered light amplitude generated at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and a signal being proportional to a scattered light amplitude generated at a second halfway point being different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way; and means for calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by a correlational function between the signal being proportional to the scattered light amplitude at the first halfway point and the signal being proportional to the scattered light amplitude at the second halfway point, the signals being observed by the means for observing the signals which are proportional to the scattered light amplitudes.
 8. The wavelength dispersion measuring apparatus for an optical waveguide according to claim 7, characterized in that the means for observing the signals being proportional to the scattered light amplitudes observes Rayleigh-scattered lights as the scattered lights at the first halfway point and the second halfway point.
 9. The wavelength dispersion measuring apparatus for an optical waveguide according to claim 7 or 8, characterized in that the means for observing the signals being proportional to the scattered light amplitudes guides scattered light generated by incident light being incident on the optical waveguide and incident light passing through delay means having delay time equal to propagation time corresponding to twice a length between the incident end of the optical waveguide and the first halfway point to first interference signal detecting means and second interference signal detecting means which have a relative delay difference equal to a distance between the first halfway point and the second halfway point to detect an interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and an interference signal I₂ being proportional to a scattered light amplitude from the second halfway point, and the means for calculating the wavelength dispersion calculates a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by a correlational function between the interference signal I₁ and the interference signal I₂ detected by means for observing the signals being proportional to the scattered light amplitudes.
 10. The wavelength dispersion measuring apparatus according to claim 9, further comprising variable delay means set on a path of the scattered light or the incident light, and characterized in that the means for observing the signals being proportional to the scattered light amplitudes respectively detects the plurality of interference signals I₁ and I₂ while causing the variable delay means to change quantity of delay, and the means for calculating the wavelength dispersion calculates a correlational function between the interference signals I₁ and I₂ by calculating the following numerical expression (6) expressing an ensemble average of products of the interference signals I₁ and I₂ detected by the means for observing the signals being proportional to the scattered light amplitudes: I₁I₂*  [Numerical Expression 6].
 11. The wavelength dispersion measuring apparatus according to claim 10, characterized in that, by using the following numerical expression (9) calculated as the correlational function of the interference signals I₁ and I₂: I₁I₂*  [Numerical Expression 8] a dispersion D in the section extending from the incident end of the optical waveguide to the arbitrary point along the way is calculated by the following numerical expression (10): $\begin{matrix} {{\frac{\overset{\_}{I_{1}I_{2}^{\star}}}{\overset{\_}{{I_{1}}^{2}}} = {\frac{\overset{\_}{I_{1}I_{2}^{\star}}}{\overset{\_}{{I_{2}}^{2}}}\mspace{59mu} = \frac{\int_{- \infty}^{\infty}{{\Gamma(Z)}{\Gamma^{\star}\left( {Z - d} \right)}{\mathbb{d}Z}}}{\int_{- \infty}^{\infty}{{{\Gamma(Z)}}^{2}{\mathbb{d}Z}}}}}{{where},{{\Gamma( z)} \cong {{\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}} z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack \begin{matrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{matrix} \right\rbrack} \right\}{\mathbb{d}\omega}}}}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 10} \right\rbrack \\ {D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}} & \; \end{matrix}$ S(ω): spectral density function of known incident light d: relative delay difference between first and second interference signal detecting means ω₀: central frequency of incident light ν_(g): group velocity of optical waveguide z₀: length of section extending from incident end of the optical waveguide to arbitrary point on the way λ: central wavelength of incident light |I₁|², |I₂|²: squares of absolute values of I₁ and I₂ β: propagation constant of output fiber.
 12. The wavelength dispersion measuring apparatus according to claim 9, further comprising variable delay means set on a path of the scattered light or the incident light, and characterized in that the means for observing the signals being proportional to the scattered light amplitudes detects the interference signals I₁ and I₂ as interference signals I while continuously change delay time τ of the variable delay means, and the means for calculating the wavelength dispersion records the interference signals I as a function I(τ) of τ and calculates a correlational function between the interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and the interference signal I₂ being proportional to a scattered light amplitude from the second halfway point by calculating the following numerical expression (8) expressing an ensemble average of the interference signals I₁ and I₂: I₁I₂*  [Numerical Expression 8] by the following numerical expression (7): I ₁ I ₂*=

I(τ_(i))I(τ_(i) +d/c)

  [Numerical Expression 7] (where <>_(i) expresses an averaging operation about τi) by using the function I(τ).
 13. A wavelength dispersion measuring method for an optical waveguide capable of measuring a wavelength dispersion in a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, comprising: a step of causing incident light having a known spectral density function S(ω) to be incident on the optical waveguide; a step of observing, by the incident light being incident on the optical waveguide in the step of causing the incident light to be incident on the optical waveguide, a signal being proportional to a scattered light amplitude generated at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and a signal being proportional to a scattered light amplitude generated at a second halfway point being different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way; a step of calculating a wavelength dispersion by a processor of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by square of correlational function between the signal being proportional to the scattered light amplitude at the first halfway point and the signal being proportional to the scattered light amplitude at the second halfway point, the signals being observed in the step of observing the signals which are proportional to the scattered light amplitudes; and a step of using the calculated wavelength dispersion to determine the characteristic of the optical waveguide.
 14. The wavelength dispersion measuring method for an optical waveguide according to claim 13, characterized in that the step of observing the signals being proportional to the scattered light amplitudes is the step of guiding scattered light generated by incident light being incident on the optical waveguide and incident light passing through delay means having delay time equal to propagation time corresponding to twice a length between the incident end of the optical waveguide and the first halfway point to first interference signal detecting means and second interference signal detecting means which have a relative delay difference equal to a distance between the first halfway point and the second halfway point to detect an interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and an interference signal I₂ being proportional to a scattered light amplitude from the second halfway point, and the step of calculating the wavelength dispersion is the step of calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by a correlational function between squares |I₁|² and |I₂|² of absolute values of the interference signals I₁ and I₂.
 15. The wavelength dispersion measuring method for an optical waveguide according to claim 14, further comprising variable delay means set on a path of the scattered light or the incident light, and characterized in that the step of observing the signals being proportional to the scattered light amplitudes is the step of detecting the interference signals I₁ and I₂ as interference signals I while continuously changing delay time τ of the variable delay means, and the step of calculating the wavelength dispersion is the step of recording the interference signals I as a function I(τ) of τ and calculating square of correlational function between the interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and the interference signal I₂ being proportional to a scattered light amplitude from the second halfway point by using the function I(τ).
 16. The wavelength dispersion measuring method for an optical waveguide according to claim 14 or 15, characterized in that the correlational function of the squares of the absolute values of the interference signals I₁ and I₂ are calculated to calculate a dispersion D in the section extending from the incident end of the optical waveguide to the arbitrary point along the way by the following numerical expression (11): $\begin{matrix} \begin{matrix} {\overset{\_}{\begin{matrix} \left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right) \\ \left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right) \end{matrix}} = {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} - \overset{\_}{{I}^{2}}}} \\ {= {K\begin{Bmatrix} {\int_{- \infty}^{\infty}{{{\Gamma(z)}} \cdot}} \\ {{{\Gamma\left( {z - d} \right)}}{\mathbb{d}z}} \end{Bmatrix}^{2}}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 11} \right\rbrack \\ {{where},} & \; \\ {{{\Gamma( z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack \begin{matrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{matrix} \right\rbrack} \right\}{\mathbb{d}\omega}}}}\mspace{11mu}\;{D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}}} & \; \end{matrix}$ S(ω): spectral density function of known incident light d: relative delay difference between first and second interference signal detecting means ω₀: central frequency of incident light ν_(g): group velocity of optical waveguide z₀: length of section extending from incident end of the optical waveguide to arbitrary point along the way λ: central wavelength of incident light K: constant |I₁|², |I₂|²: squares of absolute values of I₁ and I₂ β: propagation constant of output fiber |I|² : average of |I_(l)|² and |I₂|²  [Numerical Expression 12].
 17. A wavelength dispersion measuring apparatus for an optical waveguide capable of measuring a wavelength dispersion in a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, comprising: means for causing incident light having a known spectral density function S(ω) to be incident on the optical waveguide; means for observing, by the incident light from the means for causing the incident light on the optical waveguide, a signal being proportional to a scattered light amplitude generated at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and a signal being proportional to a scattered light amplitude generated at a second halfway point being different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way; and means for calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by square of correlational function between the signal being proportional to the scattered light amplitude at the first halfway point and the signal being proportional to the scattered light amplitude at the second halfway point, the signals being observed by the means for observing the signals which are proportional to the scattered light amplitudes.
 18. The wavelength dispersion measuring apparatus according to claim 17, characterized in that the means for observing the signals being proportional to the scattered light amplitudes guides scattered light generated by incident light being incident on the optical waveguide and incident light passing through delay means having delay time equal to propagation time corresponding to twice a length between the incident end of the optical waveguide and the first halfway point to first interference signal detecting means and second interference signal detecting means which have a relative delay difference equal to a distance between the first halfway point and the second halfway point to detect an interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and an interference signal I₂ being proportional to a scattered light amplitude from the second halfway point, and the means for calculating the wavelength dispersion calculates a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by calculating a correlational function between squares |I₁|² and |I₂|² of absolute values of the interference signals I₁ and I₂ detected by the means for observing the signals being proportional to the scattered light amplitudes.
 19. The wavelength dispersion measuring apparatus according to claim 18, further comprising variable delay means set on a path of the scattered light or the incident light, and characterized in that the means for observing the signals being proportional to the scattered light amplitudes detects the interference signals I₁ and I₂ as interference signals I while continuously changing delay time τ of the variable delay means, and the means for calculating the wavelength dispersion records the interference signals I as a function I(τ) of τ and calculates square of correlational function between the interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and the interference signal I₂ being proportional to a scattered light amplitude from the second halfway point by using the function I(τ).
 20. The wavelength dispersion measuring apparatus according to claim 18 or 19, characterized in that the means for calculating the wavelength dispersion calculates the correlational function of the squares of the absolute values of the interference signals I₁ and I₂ to calculate a dispersion D in the section extending from the incident end of the optical waveguide to the arbitrary point along the way by the following numerical expression (13): $\begin{matrix} \begin{matrix} {\overset{\_}{\begin{matrix} \left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right) \\ \left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right) \end{matrix}} = {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} - \overset{\_}{{I}^{2}}}} \\ {= {K\begin{Bmatrix} {\int_{- \infty}^{\infty}{{{\Gamma(z)}} \cdot}} \\ {{{\Gamma\left( {z - d} \right)}}{\mathbb{d}z}} \end{Bmatrix}^{2}}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 13} \right\rbrack \\ {{{where},\;{{\Gamma( z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack \begin{matrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{matrix} \right\rbrack} \right\}{\mathbb{d}\omega}}}}}{D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}}} & \; \end{matrix}$ S(ω): spectral density function of known incident light d: relative delay difference between first and second interference signal detecting means ω₀: central frequency of incident light ν_(g): group velocity of optical waveguide z₀: length of section extending from incident end of the optical waveguide to arbitrary point along the way λ: central wavelength of incident light K: constant |I₁|², |I₂|²: squares of absolute values of I₁ and I₂ β: propagation constant of output fiber |I|² : average of |I_(l)|² and |I₂|²  [Numerical Expression 14].
 21. A non-transitory computer readable medium storing a measuring program used in a wavelength dispersion measuring apparatus for an optical waveguide capable of measuring a wavelength dispersion in a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, wherein the wavelength dispersion measuring apparatus for an optical waveguide according to claim 20 is caused to receive an interference signal I₁ being proportional to a scattered light amplitude generated by incident light being incident on the optical waveguide at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and an interference signal I₂ being proportional to a scattered light amplitude generated by the incident light being incident on the optical waveguide at a second halfway point different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way as numerical values and calculate the correlational function of the squares of the absolute values of the interference signals I₁ and I2, thereby to calculate a wavelength dispersion D in the section extending from the incident end of the optical waveguide to the arbitrary point along the way by the following numerical expression (15): $\begin{matrix} \begin{matrix} {\overset{\_}{\begin{matrix} \left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right) \\ \left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right) \end{matrix}} = {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} - \overset{\_}{{I}^{2}}}} \\ {= {K\begin{Bmatrix} {\int_{- \infty}^{\infty}{{{\Gamma(z)}} \cdot}} \\ {{{\Gamma\left( {z - d} \right)}}{\mathbb{d}z}} \end{Bmatrix}^{2}}} \end{matrix} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 15} \right\rbrack \\ {{{where},\;{{\Gamma( z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack \begin{matrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{matrix} \right\rbrack} \right\}{\mathbb{d}\omega}}}}}{D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}}} & \; \end{matrix}$ S(ω): spectral density function of known incident light d: relative delay difference between first and second interference signal detecting means ω₀: central frequency of incident light ν_(g): group velocity of optical waveguide z₀: length of section extending from incident end of the optical waveguide to arbitrary point along the way λ: central wavelength of incident light K: constant |I₁|², |I₂|² squares of absolute values of I₁ and I₂ β: propagation constant of output fiber |I|² : average of |I_(l)|² and |I₂|²  [Numerical Expression 16].
 22. A non-transitory computer readable medium storing a measuring program used in a wavelength dispersion measuring apparatus for an optical waveguide capable of measuring a wavelength dispersion in a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, wherein the wavelength dispersion measuring apparatus for an optical waveguide is caused to receive an interference signal I₁ being proportional to a scattered light amplitude generated by incident light being incident on the optical waveguide at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and an interference signal I₂ being proportional to a scattered light amplitude generated by the incident light being incident on the optical waveguide at a second halfway point different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way as numerical values and calculate the correlational function of the interference signals I₁ and I₂, thereby to calculate a wavelength dispersion D in the section extending from the incident end of the optical waveguide to the arbitrary point along the way by the following numerical expression (17): $\begin{matrix} {{\frac{\overset{\_}{I_{1}I_{2}^{\star}}}{\overset{\_}{{I_{1}}^{2}}} = {\frac{\overset{\_}{I_{1}I_{2}^{\star}}}{\overset{\_}{{I_{2}}^{2}}}\mspace{56mu} = \frac{\int_{- \infty}^{\infty}{{\Gamma(Z)}{\Gamma^{\star}\left( {Z - d} \right)}{\mathbb{d}Z}}}{\int_{- \infty}^{\infty}{{{\Gamma(Z)}}^{2}{\mathbb{d}Z}}}}}{{where},{{\Gamma( z)} \cong {{\exp\left\lbrack {{\mathbb{i}}\frac{\omega_{0}}{v_{p}} z} \right\rbrack} \cdot {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack \begin{matrix} {{\frac{\beta^{''}z_{0}}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{matrix} \right\rbrack} \right\}{\mathbb{d}\omega}}}}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 17} \right\rbrack \\ {D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}} & \; \end{matrix}$ S(ω): spectral density function of known incident light d: relative delay difference between first and second interference signal detecting means ω₀: central frequency of incident light ν_(g): group velocity of optical waveguide z₀: length of section extending from incident end of the optical waveguide to arbitrary point along the way λ: central wavelength of incident light K: constant |I₁|², |I₂|² squares of absolute values of I₁ and I₂ β: propagation constant of output fiber.
 23. A method capable of measuring a wavelength dispersion of a section extending from an incident end of an optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, comprising: a step of causing incident light having a known spectral density function S(ω) to be incident on the optical waveguide; a step of observing, by the incident light being incident on the optical waveguide in the step of causing the incident light to be incident on the optical waveguide, a signal being proportional to a scattered light amplitude generated at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and a signal being proportional to a scattered light amplitude generated at a second halfway point being different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way; and a step of calculating by a processor a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by square of correlational function between the signal being proportional to the scattered light amplitude at the first halfway point and the signal being proportional to the scattered light amplitude at the second halfway point, the signals being observed in the step of observing the signals which are proportional to the scattered light amplitudes, and characterized in that the step of observing the signals being proportional to the scattered light amplitudes is the step of, when scattered light generated by incident light being incident on the optical waveguide and incident light passing through delay means having delay time equal to propagation time corresponding to twice a length between the incident end of the optical waveguide and the first halfway point are guided to first interference signal detecting means and second interference signal detecting means which have a relative delay difference equal to a distance between the first halfway point and the second halfway point to detect an interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and an nterference signal I₂ being proportional to a scattered light amplitude from the second halfway point, detecting the interference signals I₁ and I₂ through a dispersion medium which is inserted on a path on which the scattered lights from the first halfway point and the second halfway point are guided to the first interference signal detecting means and the second interference signal detecting means or a path on which the incident light is guided to the first interference signal detecting means and the second interference signal detecting means in advance and in which a value of the following numerical expression (18): β_(ref)″L_(ref)  [Numerical Expression 18] expressing an accumulated wavelength dispersion is known, and the step of calculating the wavelength dispersion is the step of calculating a correlational function between squares of absolute values of the interference signals I₁ and I₂ to calculate a dispersion D in the section extending from the incident end of the optical waveguides to the arbitrary point along the way by the following Numerical Expression (19): $\begin{matrix} {\begin{matrix} {\overset{\_}{\begin{matrix} \left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right) \\ \left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right) \end{matrix}} = {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} - \overset{\_}{{I}^{2}}}} \\ {= {K\begin{Bmatrix} {\int_{- \infty}^{\infty}{{{\Gamma(z)}} \cdot}} \\ {{{\Gamma\left( {z - d} \right)}}{\mathbb{d}z}} \end{Bmatrix}^{2}}} \end{matrix}{{where},{{\Gamma( z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack \begin{matrix} {{\frac{\left( {{\beta^{''}z_{0}} - {\beta_{ref}^{''}L_{ref}}} \right)}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{matrix} \right\rbrack} \right\}{\mathbb{d}\omega}}}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 19} \right\rbrack \\ {D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}} & \; \end{matrix}$ S(ω): spectral density function of known incident light d: relative delay difference between first and second interference signal detecting means ω₀: central frequency of incident light ν_(g): group velocity of optical waveguide z₀: length of section extending from incident end of the optical waveguide to arbitrary point along the way λ: central wavelength of incident light K: constant |I₁|², |I₂|² squares of absolute values of I₁ and I₂ β: propagation constant of output fiber |I|² : average of |I_(l)|² and |I₂|²  [Numerical Expression 20] a step of using the calculated wavelength dispersion to determine the characteristic of the optical waveguide.
 24. A wavelength dispersion measuring apparatus for an optical waveguide capable of measuring a wavelength dispersion of a section extending from an incident end of the optical waveguide to an arbitrary point along the way by using a scattering phenomenon in which a scattering coefficient at a local point on the optical waveguide is independent of time, comprising: means for causing incident light having a known spectral density function S(ω) to be incident on the optical waveguide; means for observing, by the incident light from the means for causing the incident light to be incident on the optical waveguide, a signal being proportional to a scattered light amplitude generated at a first halfway point included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way and a signal being proportional to a scattered light amplitude generated at a second halfway point being different from the first halfway point and included in the section extending from the incident end of the optical waveguide to the arbitrary point along the way; and means for calculating a wavelength dispersion of the section extending from the incident end of the optical waveguide to the arbitrary point along the way by square of correlational function between the signal being proportional to the scattered light amplitude at the first halfway point and the signal being proportional to the scattered light amplitude at the second halfway point, the signals being observed in the means of observing the signals which are proportional to the scattered light amplitudes, and characterized in that the means for observing the signals being proportional to the scattered light amplitudes detects, when scattered light generated by incident light being incident on the optical waveguide and incident light passing through delay means having delay time equal to propagation time corresponding to twice a length between the incident end of the optical waveguide and the first halfway point are guided to first interference signal detecting means and second interference signal detecting means which have a relative delay difference equal to a distance between the first halfway point and the second halfway point to detect an interference signal I₁ being proportional to a scattered light amplitude from the first halfway point and an interference signal I2 being proportional to a scattered light amplitude from the second halfway point, the interference signals I₁ and I₂ through a dispersion medium which is inserted on a path on which the scattered lights from the first halfway point and the second halfway point are guided to the first interference signal detecting means and the second interference signal detecting means or a path on which the incident light is guided to the first interference signal detecting means and the second interference signal detecting means in advance and in which a value of the following Numerical Expression (21): β_(ref)″L_(ref)  [Numerical Expression 21] expressing an accumulated wavelength dispersion is known, and the means for calculating the wavelength dispersion calculates a correlational function between squares of absolute values of the interference signals I₁ and I₂ to calculate a dispersion D in the section extending from the incident end of the optical waveguides to the arbitrary point along the way by the following Numerical Expression (22): $\begin{matrix} {\begin{matrix} {\overset{\_}{\begin{matrix} \left( {{I_{1}}^{2} - \overset{\_}{{I}^{2}}} \right) \\ \left( {{I_{2}}^{2} - \overset{\_}{{I}^{2}}} \right) \end{matrix}} = {\overset{\_}{{I_{1}}^{2}{I_{2}}^{2}} - \overset{\_}{{I}^{2}}}} \\ {= {K\begin{Bmatrix} {\int_{- \infty}^{\infty}{{{\Gamma(z)}} \cdot}} \\ {{{\Gamma\left( {z - d} \right)}}{\mathbb{d}z}} \end{Bmatrix}^{2}}} \end{matrix}{{where},{{\Gamma( z)} \cong {\int_{- \infty}^{\infty}{{S(\omega)}\exp\left\{ {{\mathbb{i}}\left\lbrack \begin{matrix} {{\frac{\left( {{\beta^{''}z_{0}} - {\beta_{ref}^{''}L_{ref}}} \right)}{2}\left( {\omega - \omega_{0}} \right)^{2}} -} \\ {\frac{\omega - \omega_{0}}{v_{g}}z} \end{matrix} \right\rbrack} \right\}{\mathbb{d}\omega}}}}}} & \left\lbrack {{Numerical}\mspace{14mu}{Expression}\mspace{14mu} 22} \right\rbrack \\ {D = {{- \frac{2\pi\; c}{\lambda^{2}}}\beta^{''}}} & \; \end{matrix}$ S(ω): spectral density function of known incident light d: relative delay difference between first and second interference signal detecting means ω₀: central frequency of incident light ν_(g): group velocity of optical waveguide z₀: length of section extending from incident end of the optical waveguide to arbitrary point along the way λ: central wavelength of incident light K: constant |I₁|², |I₂|² squares of absolute values of I₁ and I₂ β: propagation constant of output fiber |I|² : average of |I_(l)|² and |I₂|²  [Numerical Expression 23]. 